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PREFACE 


The  AGARD  Subcommittee  on  Aeroelasticity  heard  two  short  presentations  during 
the  Spring  198 1  meeting  in  £e$me,  Turkey. 

Mr  Dat  of  ONERA  presented  a  lucid  exposition  of  the  ground  vibration  testing  of 
aircraft  with  active  control  systems.  Dr  Freymann  of  Luxembourg  (and  DFVLR)  presented  a 
method  to  perform  and  interpret  dynamic  tests  on  nonlinear  systems,  such  as  control  surface 
actuators. 

The  Subcommittee  unanimously  agreed  to  publish  the  two  papers  in  one  AGARD  report 
since  they  are  an  important  contribution  to  the  technology  of  active  control. 


— 

JAMES  J.OLSEN 
Chairman 

Subcommittee  on  Aeroelasticity 
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IDENTIFICATION  DES  MODES  PROPRES  D'UNE  STRUCTURE  A  PARTIR  DES  REPONSES 
A  UNE  EXCITATION  NON  APPRQPRIEE 


par  RoIIand  DaC 

Office  National  d'Etudes 
92320  - 


et  Philippe  Dunoyer 

e£  de  Recherches  Adrospatiales  (ONERA) 
Chacillon  (France) 


REbUME 


On  prdsente  deux  mdthodes  qui  pemettent  de  traiter  lea  rdponses  harmoniques 
d'une  structure  en  un  grand  nombre  de  points  pour  determiner  les  modes  propres  de  vi¬ 
bration,  sans  que  l'excitation  ait  besoin  d'etre  appropride.  La  premifere  de  ces  md- 
thodes,  plus  synthdtique,  determine  1' ensemble  des  caractdristiques  modales  par  re¬ 
solution  d'un  systdme  d 'equations algdbriques  ;  la  seconde,  plus  analytique,  procAde  A 
une  isolation  de  chacun  des  modes  en  combinant  les  rdponses  harmoniques  des  different a 
points  de  la  structure, 

Les  rdsultats  montrent  que  ces  mdthodes  permettent  d'obtenir  une  moddlisation 
satisfaisante  de  la  structure  A  travers  la  representation  modale  , 


SUMMARY 


IDENTIFICATION  OF  THE  EIGEN  MODES  OF  A  STRUCTURE  FROM  THE  RESPONSES 
TO  A  NON  APPROPRIATED  EXCITATION 


The  two  methods  which  are  presented  use  the  harmonic  responses,  at  a  large 
number  of  points  on  the  structure,  in  order  to  determine  the  eigen  modes,  without 
necessitating  an  appropriation  of  the  excitation  forces.  The  first  one  of  the  methods, 
which  is  more  synthetic,  determines  the  set  of  modal  characteristics  through  the  re¬ 
solution  of  a  set  of  algebraic  equations  ;  the  second  one,  which  1b  more  analytic, 
proceeds  to  a  separation  of  each  mode  by  combining  the  harmonic  responses  of  the  dif¬ 
ferent  points  of  the  structure. 

The  results  show  that  the  two  methods  give  a  satisfactory  modelisation  of 
the  structure  through  ihe  modal  representation. 
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INTRODUCTION 


Les  rndthodes  d'esial  de  vibration  viaant  A  determiner  lea  modea  proprea  d'une  structure 
peuvent  ae  diviaer  en  deux  categories  :  cellea  qui  paaaent  explicitement  ou  non  par  1' appropriation  de 
l'excitation  et  cellea  qui  ne  ndcessitent  paa  1' appropriation.  A  l'dpoque  ou  l'on  ne  dlapoaait  paa  de 
ooyena  de  calcul  rapidea  aur  lea  lieux  dea  eaaaia,  il  fallait  iaoler  chacun  dea  modea  de  la  atructure 
par  une  excitation  appropride  avant  de  pouvoir  determiner  sea  caracteristiques,  [1]  A  [  3  )  .  Mais 

1'avAnement  dea  ordinateura  rapidea  et  peu  enccmbranta  offre,  depuis  quelquea  anndes ,  la  poaaibilite  de 
ddvelopper  et  d'appliquer  de  nouvelles  rndthodes  dans  leaquellea  les  rdponses  vibratoirea  de  la  structure 
subiasent  un  traitement  complexe  impliquant  le  calcul  de  systAmes  A  plusieurs  degrds  de  libertd. 

Certaines  de  ces  rndthodes  partent  de  plusieurs  eaaaia  effectuds  avec  dea  excitations  non 
appropridea  mais  inddpendantes  et  font  un  traitement  par  le  calcul  qui  revient  A  rechercher,  de  faqon 
explicite  ou  non,  la  combinaiaon  d'easais  non  approprids  simulant  au  mieux  un  essai  approprid.  [A]  A  [6], 

Ces  rndthodes  appartiennent  done  A  la  catdgorie  de  cellea  qui  paaaent  par  1' appropriation 
de  l'excitation.  En  principe  elles  ne  peuvent  dtre  appliqudea  que  ai  les  excitations  inddpendantes  atta- 
quent  toua  les  degrds  de  libertd  de  la  structure  de  maniAre  A  permettre  la  reconstitution,  par  combinai- 
son  lindaire,  d'une  excitation  appropride.  En  pratique  il  auffit  qu 'elles  attaquent  les  degrds  de  libertd 
intervenant  de  faqon  significative  dans  la  game  de  frdquence  considdrde. 

Les  rndthodes  qui  a 'affranchissent  rdellement  de  1 'appropriation  permettent  de  ddterminer 
les  caractdristiques  de  plusieurs  modes  avec  une  seule  repartition  desforces  d'excitation.  De  ce  fait 
elles  apportent  une  simplification  considerable  aur  le  plan  expdrimental  et  elles  offrent  la  possibilitd 
de  rdaliaer  des  esaais  signif icatifs  dans  les  cas  oh  la  structure  ne  peut  dtre  excitde  qu'en  un  nombre 
limitd  de  points,  voire  mdme  en  un  seul  point. 

Le  lissage  des  fonctions  de  transfert  par  des  fractions  rationnelles ,  devenu  classique 
depuis  quelques  anndea,  eat  une  de  ces  rndthodes.  11  ne  sera  pas  ddcrit  ici,  parce  qu'il  a  ddJA  fait  l'ob- 
jet  de  plusieurs  publications,  17]  A  [10]  ,  et  nous  nous  bomerons  A  rappeler  qu'A  chaque  degrd  de 

libertd  (ou  A  chaque  point)  de  la  structure,  il  fait  correspondre  une  fraction  rationnelle  dont  les  pSles 
s ' identif ient  aux  pOles  de  la  atructure  elle-mdme,  tandis  que  lea  rdsidus  ddfinissent  les  formes  proprea. 

Le  lissage  par  fraction  rationnelle  a  dtd  appliqud  par  l'ONERA  A  une  structure  du  lanceur 
ARIANE  excitde  en  une  seule  section.  Les  modes  sont  ddterminds  de  faqon  satisfaiaante,  du  moina  pour  les 
points  de  grande  amplitude,  mais  les  temps  de  calcul  sont  assez  longs  et  la  mdthode  prdsente  l'inconvd- 
nient  de  donner  des  pOles  entachda  d'une  erreur  expdrimentale  diffdrente  d'un  point  A  1 'autre  de  la  struc¬ 
ture. 


Le  prdsent  exposd  est  limitd  A  deux  rndthodes  qui  ne  prdsentent  pas  ce  dernier  inconvdnient , 
car  elles  prennent  en  considdration  1 'ensemble  des  points  de  la  structure  pour  ddterminer  chacun  des  pOles 
et  ne  fouraissent  done  qu'un  seul  pOle  par  mode.  La  premiAre  de  ces  rndthodes,  plus  synthdtique,  ddtermine 
l'ensemble  des  caractdristiques  modales  par  rdsolution  d'un  systAme  d'dquations  algdbriques  ;  la  seconde, 
plus  analytique,  procAde  A  une  isolation  de  chacun  des  modes  en  combinant  les  rdponses  harmoniques  des 
diffdrents  points  de  la  structure. 

Ces  rndthodes  peuvent  s'appliquer  A  n'importe  quel  systAme  rdgi  par  un  systAme  d'dquations 
dif fdrentielles  lindaires  A  coefficients  constants,  mais  e'est  uniquement  1 'application  aux  structures 
qui  sera  ddveloppde  et  illustrde  ici. 

Remarquons,  pour  terminer  cette  introduction,  qud  si  ces  rndthodes  s ' acconmodent  d'une  exci¬ 
tation  non-appropride,  il  ne  faut  pas  en  ddduire  que  l'excitation  peut  dtre  quelconque.  La  ddtermination 
d'un  mode  n'est  possible,  en  effet,  que  si  l'excitation  lui  fournit  assez  d'dnergie  pour  qu'il  dmerge 
suf fisament  dea  rdponses  vibratoirea.  A  ddfaut  d'dtre  appropride,  l'excitation  doit  done  posseder  un 
minimum  d'efficacitd  pour  les  modes  que  l'on  ddaire  identifier. 

I  -  METHODE  DE  SYN THESE  MODALE 


La  structure  est  assimilde  A  un  systAme  lindaire  discret  A  tv  degrds  de  libertd  rdgi  par 
l'dquation  du  second  ordre  : 


MX+BX+KX  =  t 

I-I-I 

ou  M  ,  B  et  K  sont  des  matrices  rdelles,  symdtriques,  positives,  A  coefficients  constants, 
de  dimension  x  n. 

X  et  E  sont  des  colonnes  de  dimension  Tl  qui  ddpendent  du  temps  :  X(l)  est  le  vecteur 
rdponse  d»  systAme  et  le  vecteur  excitation. 

Cn  admettra  que  les  variables  X  qui  constituent  le  vecteur  X  sont  accessibles  A  la 
mesure  et  qu'elles  sont  inddpendantes.  "  t  t 

En  rdgime  harmonlque,  X(k)  et  EftJ  sont  remplacda  par  g  et  E  0  ,  X  et 

C  ddsignant  dans  ce  caa  dee  colonnea  inddpendantes  de  t.  On  a  : 


[.  +  ao,  B  +  K  ]X  ~  E. 


I-I-2 


I 


3 


La  co tonne  E  e«t  rielle  si  1' excitation  est  AqulphasAe  et  prut  pour  rifirence  de  phase. 
Dan*  ce  qui  *uit  on  admettra  qua  cette  condition  eit  rialisie  dan*  le*  e*tai*  considiris. 

Au  court  de  l'etaai  hanaonique,  l'excitation  £  eat  imposie  et  l'on  determine  le  vecteur  riponse  X 
A  partir  dea  nature*  de  diflexion.  On  conaidhre  lea  rdaultata  obtenua  lore  d'un  balayage  en  frequence  ef- 

£  indipendant  de  la  pulaation  oJ  . 

.i 

M  K 


fectui  avec  un  vecteur  excitation 

I -I  Ditemlnation  de*  produita  natriclela 

Compte  tenu  de  la  aymitrie  dea  matrices  M  , 


MB 


s 'Acrirt 


et 

«  ic 


.a 

M  £ 


1' Aquation  (I-I-2)  peut 


KM 


EH'1  sont  inconnues 


[X  -i] 

Dans  cette  Aquation,  le  vecteur  X  (iM))  eat  connu  et  lea  matrices  KM*" 


BM* 

EM-1 


2  rz 

-  co  X 


1-1-3 


BM''  et 


L'Aquation  (I-I-3)  ne  peut  pas  Atre  rAsolue  directement.  De  plus,  dans  la  plupart  des  ap¬ 
plications  pratiques,  elle  ne  peut  Atre  virifiie  que  de  faqon  approchAe  A  toutes  les  valeurs  de  Cj  , 
parce  que  la  mesure  du  vecteur  X  est  entachAe  d'erreurs  et  que  la  structure  essayie  n'est  pas  rigou- 
reusement  confome  au  noddle  liniaire,  discret,  du  second  ordre. 

Une  solution  approchAe  peut  Atre  obtenue,  pour  chacune  des  colonnes  de  la  matrice  inconnue, 
en  cherchant  le  minimum  du  scalaire  £r  dAfini  par  : 


J[x  :-«-■*  ’-•] 

-00 


KM* 

BM'1 

lEK1] 


e  .  cfx  e 

r 


fr  Id£° 


OU  c„  est  le  vecteur  cl4  qui  extralt  la  r*  colonne.  En  r€solvant  ce  probl&me  de  minimum  comae  indi- 
qu€  dans  1' annexe  A,  pour  toutes  les  colonnes  (  t  *1,2, . . . TL  ) ,  on  obtient  la  solution  : 


0 0 


KM'1 

/D  f 

X 

BM1 

= 

(K 

AOJX 

.E  M‘\ 

j 

0 

-1 

[x\-uoX\ -i] 


dio 


X 

iu>X 

-1 


X  dui 


L'Aquation  (1-2)  ditermine  done  la  matrice  inconnue  A  partir  du  vecteur 
duit  les  produits  matriciels  M  K  ,  M*  B  cr  ME  en  transposant. 


1-2 

(4M).  On  dA- 


En  pratique,  les  intigrales  de  0  A  l'infini  peuvent  Atre  remplacAes  par  des  sosanations  sur 
les  frAquences  de  mesure. 

1-2  DAtermlnation  des  caractArlstlques  modales. 

La  reprAsentation  modale  eBt  rialisie  en  effectuant  le  changement  de  variable  qui  rend 
les  matrices  de  masse  et  de  raideur  diagonalea.  Ce  changement  de  variable  est  dAfini  par  une  matrice  in¬ 
connue  Q  ,  carrAe,  riguliire,  de  dimension  *V1  >  T1  .  On  a  : 

x  =  0a 

est  la  colonne  des  variables  de  la  reprAsentation  modale. 


1-3-1 


oh 


En  subatituant  dans  (1-2)  et  en  primultipliant  par  @ 

[-wyU  ^  =  Q 

pSM0 

y=SKe 

©B© 


il  vient: 


1-3-2 

1-3-3 

1-3-4 

1-3-5 


11  faut  diterminer  la  matrice  0  qui  rend  JA  et  ^  diagonales,  apris  quoi  on  devra 
diterminer  les  matrices  yU  et  £  ellea-mAmes  et  la  matrice  ^  . 

Pour  ditermlner  0  ,  on  part  du  produit  matriciel  M  K  diduit  de  l'Aquation  (1-2). 


1  _  .1 

D’aprds  (1-3-3)  et  (1-3-4)  on  a:  M"  K  =  &  f*  Q  & 

.1 

Le  produit  matriciel  u  V  eat  la  dlagonale  dea  valeurs  proprea  du  ays t feme  conservatif  associd.  En  poaant 

19  2  -1  - 

/  1-4-1 

on  ddduit  que  lea  pulaatlona  proprea  et  lea  vecteura  ©e.  de  la  matrlce  0  aont  lea  solutions  du 
probldme  :  *  * 

[  v2-  m  'k  J  Qe  =  o 

r  r  1-4-2 

Connaissant  la  matrice  et  le  produit  matriciel  mV  ddduit  de  l'dquation  (1-2),  on  determine  le 

produit  matrlclel  d'aprfes  (1-3-3)  et  (1-3-5)  : 

yf'fl  =  ©‘MB  © 

Enfin,  pour  determiner  la  matrlce  jA  ,  il  faut  partlr  des  vecteurs  mV  ddduits  de  l'dquation  (1-2).  On 

mV  - 

d’°*  :  H  mV  =  £  E 

d'oh,  dtant  dlagonale,  et  en  ddsignant  par  l'dldment  Q  jj(  £ 


u  - 

A  er  s'ff't 

er©E  ddaigne  le  r'  ^  terme  du  produit  matriciel  0£  ,  et  de  m&ae  pour  ereW’E. 

t  On  remarque  que  la  matrice  de  changement  de  variablea  e  et  les  produits  matriciels  if  ^ 

et  /*(>  aont  ddterminds  4  partlr  du  vecteur  rdponse,  X(ioS)  ,  sans  qu'il  aoit  n£cessaire  de  connaiti 


et  aont  ddterminds  4  partlr  du  vecteur  rdponse,  X(ico)  ,  sans  qu'il  aoit  ndcessalre  de  connaitre 

l'excitation  fc  .  Celle-ci  n'intervient  que  dans  la  formule  (1-4-4)  pour  le  calcul  des  masses  gdndrali- 
ades  it 

A 

2  -  METHP  .E  D' ISOLATION  MODALE 

2-X  Caa  gdndral.  Dans  la  mdthode  prdsentde  dans  ce  paragraphe,  1 'Isolation  des  modes  est  rdalisde  en 
auperpoaant  lindairement  lea  rdponses  des  diffdrents  degrds  de  liberty  4  une  excitation  harmonique. 

Comma  dans  le  paragraphe  prdcddent,  la  structure  eat  assimilde  4  un  systfeme  lindaire  diacret 
4  n.  degrds  de  libertd  rdgi  par  l'dquation  : 

MX  +  BX  +  KX =  £ 


On  obtlent  la  forme  canonique  en  posant  : 


II  vient  : 


(i«x) 

■^-[.rr'x  -m-'bJ 

r-  (A) 

Z  (iw)  =  [iw  -  A]  f 


La  matrice  A  eat  carrde,  de  dimension  !  n  x  2H  ; 
Z  et  aont  des  colonnea  de  dimension  2  TV  ; 

Z  (tfll)  eat  le  vecteur  d'dtat  du  syatfeme. 

2-1-1  Principe  de  la  adparation  modale 


Suppoaona  que  .A  posadde  2  valeurs  proprea  diatinctea  y  ,  et  ddaignona  par'll  et  lea  vecteura 

proprea,  respectlvement  4  drolte  et  4  gauche,  de  A  .  Lea ^vecteura  K  et  Hi  aont  orthogonsux,  pour 
l  It  ,  et  ai,  de  plus,  lea  vecteura  aont  convenablement  normalisds,  ou  a  : 


V  u  =  X 

w  i  wi 


1 


(  symbole  de  Kronecker) 


5 


Lei  «  ,  U,  et  V  constituent  lea  d laments  de  la  decomposition  en  elements  simples  du  vecteur  z  (*<*) 


k  XU-, 


2-6 


Ed  prdmultiplisnt  par  l'un  des  vecteurs  Vt  ,  pula  en  transpoaant,  on  obtient  I'dquation  : 

2TM  V  =  -£.?£- 

*’ui~  Tr 

Cette  Equation  montre  que  la  superposition  des  variables  d'dtat  ddfinie  par  le  vecteur  1JZ  isole  un  mode 
propre  du  systdme. 

Connaiasant  le  vecteur  d’dtat  2L  (qui  peut  6tre  ddduit  des  mesures  de  rdponse  dans  un  essai  hannonique 
avec  balayage  en  frequence),  on  peut  trouver  les  solutions  de  (2-6),  c'est-&-dire  les  et  les  ■*»  , 

par  iteration.  r 

2-1-2  Recherche  des  solutions  (>U)  par  iteration 
Les  vecteurs  %  peuvent  6tre  normalises  de  faqon  arbitraire.  Si  on  les  normalise  en  posant  : 


1 'equation  (2-6)  devient  : 


y  %=  i 

zT(i*y vr  = 


2-7 


2-8 


A.tO-fr 

Les  solutions  de  cette  equation  peuvent  s'obtenir  par  iteration  :  partant  d'une  valeur  estimee  de  jl»  on 
determine  le  vecteur  “\fr  qui  vdrifie  I'equation  (2-8)  de  faqon  exacte  ou  approehde,  dans  une  ganine  de 
frequence  assez  large,  puis  partant  de  ce  vecteur  on  determine  ^  et  on  recommence. 

Mais  lorsque  et  sont  complexes,  il  est  plus  avantageux  de  combiner  liquation  (2-8) 

**"  d^signant  par  l1  indice  qui  rep&re  la  valeur  propre 

normalisation  (2-7)  et  1  Equation  ^2-8)  devien- 


relative  k  ^  avec  1 'Equation  relative  k  J*  .  En  design 
•p*  et  le  vecteur  propre  correspondent,  ti  condition  de 
n4nt  :  7- 

7  IT,  =  ' 

zhv;  - 


4.0). 


X 


D’ou,  en  faisant  la  somme  et  la  difference  avec  (2-8)  : 


En  posant  : 


il  vient  : 


ZThr  +  IT]  =  l^-CTr^rr) 


zTrv  vi-  xcJt 

l  r-  -r J  .o>2_W^ 


+  t r  )  +  TrTr 

fr=  /Vr- *rVr 

t  =-#  Cif-vr) 

TT  ^  y 

^  (***)  }jr  ~  -  cu2+  2-CWc^V  +  v/0  +  °^) 


2-9-1 


2-9-2 


2-9-3 


2-10-1 


^  +  ^rVr 


-  (j/+  2it0  of.  Vf  +  V  0  +  °<t) 

On  peut  determiner  les  solutions  de  (2-I0-I),  c'eat-4-dire  IT  o(  et  V.  ,  par 
0/  Vf  dans  (2-10-2)  et  celculer  ^  .  #f  f  1 


2-10-2 


tdratlon,  puis  reporter 


Pour  les  raisons  d<)S  invoqudes  (erreurs  de  mesure  et/ou  systdme  non  rigoureusement  conforme  au  noddle 
thdorique),  les  equation*  (2-8)  et  (2-10)  n'admettent  souvent  que  des  solutions  approchdes. 


6 


Celles-ci  doivent  Cere  calculdes  &  l'aide  d' algor ithmes  de  resolution  par  lea  moindres  carrds-  Par  exemple, 
la  recherche  dea  solutions  de  (2-I0-I)  par  iteration  comporte,  S  chaque  pas  d'iteration,  deux  sequences 
(e  lissage  par  les  moindres  carrds  : 

f  -  dans  la  preaidre  sequence,  partant  des  valeurs  d  et  V  du  pas  precedent,  on  calcule  le 
vecteur  qui  rend  minimum  l'dcart  ddfini  par  : 

£  _  ^ I  2V - - - - r-fotco 

'  J  I  "t  -  lOZ  +  eitJoi 'V '  +  vVf+oQ  I 

_oo  r  r  f  v  * 

et  M  est  obtenu  en  appliquant  la  formule  Cl-1  de  1' annexe  A. 

w  $L  /  t 

-  dans  la  deuxifcme  sequence,  on  po9e  £o(rVr  =  b  }  a-  C  ,  on  introduit  une 

inconnue  auxiliaire,  Cl  ,  et  on  cherche  let'  valeurs  de  cl  fa  ei“  C  qui  rendent  minimum  £  ddfini 

nor  •  an  A  2 


-  f  |  +•  c  )Z^"~  a- 


A  chaque  paa  d1 iteration,  1' inconnue  auxiliaire  CL  doit  tendre  vers 

En  pratique,  les  integrates  de  -00  Q  Oo  peuvent  Stre  remplacdes  par  des  sommations  sur  CO  et  «-  CO  , 
les  CO,  6tant  les  frequences  de  mesure.  ^  * 

1C 

Remarque  I  -  * 


TTr  etant  le  vecteur  propre  associd  h  la  valeur  propre  ,  son  cjpnjuguS  est  vecteur  propre  as- 

socid  £  =  V_r  ♦  Le  vecteur  V.f  se  ddduit  done  du  vecteur  en  lui  appliquant  la  condition  de 

normalisation  -  1  *  e'est-fc-dire  que  l*on  a  :  ** 

-r  wv 

ryf  r 

Si  le  vecteur  Jr  est  rdel  (cas  de  l'excitation  6quiphasde  et  prise  cotnme  rdffirence  de  phase),  on  a 
?V=  1  (d'aprbs  (2-7),  d'oO  :  lT  =  IT* 

f  \  r  ’  /  -rr  j 

On  en  ddduit,  d'aprds  (2-9-2)  et  (2-9-3),  que  U  eat  la  partie  rdelle  de  V  et  II  la  partie 
imaginaire  :  *  &C 

t-  i 

fr-iV-tf) 

/  * 

Les  vecteurs  g  et  des  Equations  (2-I0-I)  et  (2-10-2)  sont  done  r€els  lorsque  J*  est  r€ei. 


Remarque  2  - 

On  voit  que  l'isolation  modale  par  iteration  peut  fitre  effectu^e  sans  qu'il  soit  n6cessaire  de  connaltre 
le  vecteur  excitation  et  la  matrice  ht 

Remarque  3  - 

Le  processus  d' isolation  modale  donnant  les  valeurs  propres  et  les  vecteurs  propres  &  gauche  % J"  , 
on  peut,  en  principe,  reconstituer  la  matrice  et  le  vecteur  . 

A  =  (u7>  vT 

0\J  IT  est  la  matrice  des  vecteurs  Vr  et  P  la  matrice  diagonale  des  valeurs  propres  ^>f  .  Le  vecteur 
«JT  est  ddtermind  par  la  condition  de  normalisation,  1JT  T*  -  I,  d’oCi  l’on  d£duit  : 

*  -I  /  4  \  r 

*-Cirr(i) 

2-2  Systfeme  passlf  faiblement  disaipatif 

La  plupart  dea  structures  peuvent  Stre  assimildes  &  des  systdmes  passifs  faiblement  dissipatifs.  Avec  la 
representation  modale,  l'dquation  de  la  dynamique  devient  quasi  diagonale  et  s'dcrlt  : 

[-"V  *iu)P +  *]  v  ^ 

oO  yU  et  y  sont  des  matrices  diagonales  et  une  matrice  symdtrique  mais  gdn^ralement  pas  diagonale. 
La  matrice  entre  crochets,  peut  6tre  cons iddr4e* comae  quasi  diagonale  parce  que  les  coefficients  /$ 
sont  petlts.  1  tS 


La  forme  cancnique  est  obtenue  en  posant  : 
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f-/V  -/ji  ] 


2-II-3 


2-1 1-4 


L'hypothfcse  du  systfeme  faiblement  dissipatif  peut  ae  traduire,  sur  le  plan  raatlidmat ique  en  conaiddrant  lea 
coefficients  de  la  matrice  yj  corane  dea  infiniment  petit8  du  premier  ordre,  alora  que  lea  coefficients  de 
II  et  de  y  sont  des  grandeurs  finies.  Si  l'on  admet,  en  outre,  que  les  frequences  propres  aont  auffiaam- 
ment  distinctea  pour  que  ^  ^  aoit  une  grandeur  finie  »  on  peut  expliciter  y  juaqu'au 

premier  ordre.  11  vient  (voir  annexe  B)  : 


/  2T  [  l  +  *  V]  e. 


2-12-1 


ou  eat  une  matrice  de  dimension  H  a*  dont  les  coefficients,  infiniment  petits,  aont  ddfinis  par  : 


4Ve.'  = 


I _ Vr  ft  r- 

Al'C-'f) 

O 


It  r  t 


2-12-2 


En  rempiaqant  J*  par  (2-II-4),  on  a  : 


*yTor  =_fr  e^(i**yT)Q 


et  le  vecteur  if.  normalise  par  la  condition  *¥  if s  1  eat  donn6 


par  : 


/  -y  U  + 

rreT[i^yTJQ 
l4-xq/]er 
\  [JL  +  *yTJQ 


Ccnaiddrons  le  cas  de  l 'excitation  dquiphas£e,  ce  qui  rend  la  colonne 
formule  (2-10-1)  au*  coordonndes  gdndralisdea  CJ  on  a  : 


r  T  .  T  -j  Mr 

M  !  A£0(i  c/r*"  -coz+2icoo(rV  +  \>.1 


2-13 

r6elle.  En  appliquant  la 

2-14 


(  a  6t6  nfigligfi,  gtant  infiniment  petit  dans  1'hypothSse  du  systSoe  faiblement  dissipatif). 

Dans  cette  Equation,  est  la  partie  imaginalre  du  vecteur  Vr  d6fini  par  (2-13).  En 

remplaqant  ■»  par  XV  _  o(  V  et  en  nigligeant  les  termes  du  second  ordre,  la  formule  (2-13)  donne  : 

r  r  r  r  t  r'  \ 

X-  c; 

r'  r* 

ou  v  et  C  sont  des  colonnes  de  dimension  fl  donndes  par  : 

C'  . 

c  v,  e,  Q 

O'.  y“  e.err- e,e,VT)  (} 
r "  (  <Q)‘ 

Maia  en  g^ndral  on  ne  peut  pas  accdder  directement  aux  coordonndea  modales  Q  et  on  meaure  d'autres  va¬ 
riables  X  aasimilables  I  des  degrds  de  libertd  :  les  vecteurs  X  et  Q  sonx  relics  par  une  Equation  de 
changement  de  variables  ddfinie  par  une  matrice  carrde  rdelle  Q 

X  =  9  c, 

Des  Equations  (2-14),  (2-15)  et  (2-16)  on  d6duit  : 


2-I5-I 


2-I5-I 


2-15-3 


2-16 


2-I7-I 


8 


oik 


f  et  Pt  sont  del  vecceun  riels  d'ordre  n.  difinii  par 

Jr 


%  -  f  # 
-  -*  , 


2-17-2 


2-17-3 


„  f0I>  Q 

V — (<q7 

Le«  coefficients  o(r  ,  >{*  et  lee  vecteure  F,  et  O.  peuvent  £tre  ddtermind*  en  resolvent  liquation 
(2*17-1)  par  iteration  canme  indiqud  dans  le  parfgraphe  (2-1-2). 

On  remarque  que  le  vecteur  f)  est  nul  si  la  matrice  de  dissipation  eat  diagonale,  c’est-d-dire 
sil'ona  /3  ;o  ,  V#*  '  r 

Determination  del  matrlcei  te  et  /M 

*  —  /y 

Lea  masses  generalities  ^Uret  la  matrice  @  se  diduisent  del  vecteurs  Jr  .  On  remarque  que  f  Sr  Q 

eat  un  noobre  riel  qui  depend  de  la  nomaliaation  des  formes  propres  et  qui  a  la  dimenaion  d'un  Cemps. 
Fixer  arbitralrement  la  valeur  de  ce  rapport  revient  done  &  noroaliaer  la  forme  propre.  Poaona  par  exeople 

AVr  J_ 

<£  0  2-18 

En  substituent  dans  (2-17-2)  il  vient  : 

f©7'er=  U 

d'ou,  en  disignant  pare  la  matrice  des  vecteurs  JT  et  par  V  la  matrice  diagonale  dea  frequences 
propres  Vf  ,  '  r 


2-19 


e  =  [jrT  W* 

Lea  masses  gindralisdes  sont  diduites  de  la  condition  de  normalisation  C 2— IB  >  : 

La  colonne  C)  eat  ddterminde  si  l'on  connate  le  vecteur  excitation  EL  dans  le  systbme  de  variables 
X  .  Le  thiorerae  des  travaux  virtuels  donne  en  effet  : 

SX  E.  —  Q  d'ou  l'on  diduit,  d'apris  (2-16)  : 

Q--GTt. 


Determination  des  coefficients 


Leg  coefficients  non  diagonaux  de  la  matrjee  fi  peuvent  Stre  diduits  des  vecteurs  >J  donnis  par  (2-17- 
En  premultipliant  ^  par  la  matrice  Q~on  a  :  'r 

&\  =  ti 
,r  «Q)‘ 

*  T 


Go, 

terme.  1 


eat  un  vecteur  d'ordre  n  .  Multiplions  le  &  gauche  par  (> 
term^.  11  vient  :  » 


■-P-  r 


pour  extraire 


I  14 
le  I r 


e:eT7  =  e<Viv/- <Q 


«  or 

A  eJVe, 

<Q 

.  —-L@i is. _ 

ttVK  <? 


(  d'aprea  2-12-2  ) 


D'oa 


A,- 


2-21 


1 
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r/t  Jr  AT  .  itim  »T 

oik  CrW  eat  le  f  coefficient  de  la  colonne  M  et  C7  le  k  coefficient  de  la  colonne 

3  -  APPLICATION  A  L'lDEHTIFICATIOH  DES  CARACTERISTIQUES  HODALES  D'UNE  STRUCTURE. 

Lea  mCthodes  prdsentdes  aux  paragraphea  I  et  2  penaettent  de  determiner  lea  caractdriatlquea 
mode lea  d'une  atructure  A  partir  d’un  eaaai  harmonique  avec  balayage  en  frequence  a  force  conatante. 

On  admettra  que  lea  meaures  de  rdponse  Bont  effectudes  en  dea  pointa  suffiaamnent  nombreux  et 
bien  rdpartia  pour  que  lea  vecteura  ddfiniaeant  la  reatriction  A  cea  pointa  dea  tv  formea  proprea  in* 
tervenant  dana  la  game  de  frequence  exp  lor  de  aoient  inddpendanta .  Cette  condition  impllque  que  toua  lea 
degrda  de  libertd  intervenant  dana  cette  game  de  frequence  aoient  acceaaiblea  A  la  meeure. 

3-1  Reduction  du  ay at feme 

Soit  H  le  nombre  de  pointa  de  meeure  et  N  le  nombre  de  pointa  d'excitation.  En  auppoaant  que  lea 
pointa  d'excitation  aont  dgalement  pointa  de  meeure,  on  a  M  .  En  pratique,  H  eat  infdrieur  ou  dgal 

au  nombre  de  degrda  de  libertd,  tv  ,  et  M  eat  aouvent  beaucoup  plua  grand  que  qv  • 

Soit  v/*(  4-Uj]  la  colonne  de  dimension  n  qui  deflate  les  deflexions  aux  points  de  mesure  et 

F  la  colonne  d'ordre  H  qui  d^finit  la  repartition  des  forces  d'excitation. 

Les  deflexions  aux  points  d'excitation  sont  d^finies  par  une  colonne  ,  de  dimension  if  , 

d^duite  de  au  moyen  d'une  matrice  de  permutation  %J 

%  *  J  Kt  3-1 

Si,  par  exemple,  il  y  a  deux  excitateura,  aux  points  d'indlce  2  et  3,  on  a  : 


j  *  t : : :  i  ■  (§) 


Base  modale,  base  artificielle 

La  restriction  de  chacune  des  formea  proprea  aux  points  de  mesure  eat  ddfinie  par  un  vecteur  <P.  ,  de 
dimension  l~I  ,  et  1* ensemble  de  cea  vecteura  forme  une  matrice  <P  de  dimension  M*n  .On  pose  : 

<p  <J  (<*>)  3-2 

eat  la  baae  modale  qu'il  faut  determiner. 

Cons iddr ant ,  d'autre  part,  tv  frequences  voiainea  dea  frequences  de  resonance  ,  on  peut 

meaurer,  pour  chacune  de  cea  frequences,  le  vecteur  5,  <Jui  dSfinit  la  composante  de  la  rdponse  en  qua¬ 
drature  par  rapport  A  l'excitation  aupposde  equiphaade. 

Cea  vecteura,  de  dimension  M  ,  torment  une  matrice  rdelle  S  ,  de  dimension  Mart  ,  qui  peut 
conatituer  une  baae  de  representation  des  deflexions,  au  meme  titre  que  .  On  a  : 

%  (<■*)  =  SX(iu))  3.3 

oil  X  est  une  colonne  de  mfene  dimension,  ^  ,  que  la  colonne 

X  eat  d6 termini  par  la  m4thode  des  moindres  carr€s  &  partir* des  meaures  de  deflexion 


X  ( tu»)  =  [5tS]Vv/(a«j) 


Lea  oethodea  dea  paragraphea  X  et  2  peuvent  fitre  appliques  au  vecteur 


Changement  de  baae. 


©  etant  la  matrice  de  changement  de  base  reliant  X  et  Cj  par  l’dquation  (1-3-1),  X-  0q 


on  dddult,  d'aprAa  (3-2)  et  (3-3) 


<f>  =  5© 

Cette  equation  tervira  A  determiner  ^  lorsqu'on  connattra  . 


Equation  de  la  dynamlque 


L'equatlon  de  la  dynamique  peut  8tre  formuiee  par  rapport  aux  variables  du  vecteur  X  ou  par  rapport  aux 
coordonndes  gendraliseea  de  la  representation  modale,  vecteur  Cj  • 

Dans  le  premier  cas  on  a  : 


+  K  ]x  =  e 

E  =  ST  Jr  F 


3-6-2 
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ec  dans  le  cas  de  la  representation  modale  : 


avec 


Q  =  0T  SrJ 


0 


3-6-3 

3-6-4 


3-2  Mdthode  de  synthfese  modale  du  paragraphe  I 


formule  (3-4),  on  peut  appliquer  la  formule 


Le  vecteur  XM  dtant  ddduit  des  mesures  de  deflexion  par  la  fornule  ,  _  _ 

(1-2)  et  determiner  les  produits  matriclels  M*' K  ,  M'1  B  et  M"1  E.  .  On  calcule  ensuite  Q  et  lea 
prodults  matriclels  yU  p  et  ,  par  les  formules  (I-4-I)  4  (1-4-3),  puls  les  masses  gdndrallsdes  JUy 

par  la  formule  (X-4-4).  Enfln  les  formes  propres  de  la  structure  qui  constituent  la  base  modale  sont  don- 
ndes  par  la  formule  (3-5). 

3-3  Isolation  des  modes,  mdthode  du  paragraphe  2 

Le  vecteur  X  u«)  etant  connu ,  on  cherche  les  solutions  fr  ,  nr 
( 2—17  —I )  en  proc^dant  par  itdration  conme  indiqud  au  paragraphe  (2-1-2). 


°< r 


la  formule  (2-19) 


,  de  l'dquation 

La  matrice  0  est  dbterminde,  4  partlr  des  vecteurs  €  et  des  pulsations  propres  Vr  ,  par 
e  (2-T9)  :  Jr 

e  -  t 


Les  masses  gdndralisdes  sont  donndes  par  les  formules  (2-20)  et  (3-6-4) 


A 


eTQTSTJTF 


3-7 


et  les  formes  propres  sont  d6termin6es  par  (3-5) 


<p  ^  SQ 

es  yK*.  et  les  ainsl  trouvds  sont  ddterminds  par  la  condition  de  normalisation  (2-18).  11 

,  apres  coup  de  se  ramener  A  une  normalisation  plus  conventionnelle. 


Les 

est  facile 

3-4  Excitation  appropride  A  un  mode  du  systfeme  conservatif  associd. 

L'excitation  est  appropride  4  un  mode  r  du  systdme  conservatif  associd  si  la  colonne  Q  donnde  par 
(3-6-4)  est  proportionnelle  au  vecteur  Br  •  Le  vecteur  F"  qui  ddfinit  la  rdpartition  de  force  appro- 
pride  dolt  done  dtre  une  solution  approchde,  ou  dventuellement  exacte,  de  l'dquation  : 

©Vj>  =  erc 

ou  C  est  une  constante  qui  a  la  dimension  d'un  travail. 

Le  nombre  de  lignes  du  produit  matriciel  6TsrJT  est  dgal  &  Yl  (nombre  de  degr€s  de  liberty) 
et  le  nombre  de  colonnes  &  M  (nombre  de  points  d 'excitation) .  Si  N  <  tv  on  a  la  solution  approchde  : 


f=  [  jsg  eTsrfJ  G  sJe^c 


et  al  N  s  H  on  a  la  solution  exacte 


f  =  ( ©T sr J  )  erc 


.1 


3-8-2 


3-8-3 


Pour  ddterminer  le  coefficient  C,  11  faut  se  donner  une  condition  portant,  par  exemple,  sur 
1 'amplitude  de  la  rdponse. 

L'appropriation  de  l'excitation  ne  peut  etre  rdallsde  que  lorsqu'on  a  ddjA  proeddd  4  une  premidre 

Identification  modale  avec  une  excitation  arbltralre.  Sans  dtre  Indispensable,  elle  permet  de  refaire 

dansde  mellleures  conditions  1 ' identification  des  modes  pour  lesquels  l'excitation  arbitraire  initiate 
pourralt  se  rdvbler  trbs  mal  adaptde. 

3-5  Rdsultats 

Les  figures  (1)  4  (3)  et  le  tableau  (1)  illustrent  les  rdsultats  donnds  par  les  mdthodes  des  psragraphes 
(I)  et  (2)  sur  l'avion  T  B  20,  monomoteur  4  aile  drolte  de  la  Socidtd  SOCATA.  Le  nombre  de  points  de  me- 
sure  prls  en  considdratlon  dans  cette  application  btait  dgal  4  20  (  M  ■  20). 

On  sxmtre  id.  les  rdsultats  obtenus  4  partlr  d'un  balayage  en  frdquence  avec  excitation  en  haut 
de  la  ddrive,  dans  la  gamme  de  10  4  20  Mz  . 

L'allure  des  courbes  de  rdsonance  montre  qu'il  exlste  quatre  modes  dans  cette  ganme  de  frdquence 

et  qu'il  faut  done  introduire  quatre  degrds  de  llbertd  dans  les  calculs  (  -tv  -4). 


Les  caractdristiques  modalea  one  dtd  ddtermindes  par  lea  mdthodes  dea  paragraphea  (I)  et  (2). 

Le  tableau  (1)  permet  de  comparer  les  frequences,  lea  amort  Easements  et  lea  masaes  gdndralisdes 
donnds  par  lea  deux  mdthodes. 


TABLEAU  I 


FREQUENCE  (HZ) 

AMORTI S  S  EMHHT 

MASSE  GENERALISES  (m2  kg) 

Mdthode  1 

Mdthode  2 

Mgthode  1 

Mdthode  2 

Mdthode  1 

Mdthode  2 

Mode 

1 

11,76 

11,77 

0,020 

0,028 

25,62 

16,99 

Mode 

2 

14,31 

14,33 

0,029 

0,030 

17,91 

15,09 

Mode 

3 

15,59 

15,59 

0,008 

0,012 

172,7 

152,4 

Mode 

4 

17,49 

17,49 

0,038 

0,041 

12,56 

10,93 

La  figure  (2)  montre  la  comparaison  dea  courbes  de  resonance  expdrimentales  et  des  courbes 
reconstltuees  par  le  calcul  A  partlr  dea  caracterlatlques  modalea  donnees  par  lea  deux  methodes.  La  posi¬ 
tion  aur  l'avion  dea  points  ou  sont  effectuds  ces  comparaiaons  eat  ddfinie  dans  la  figure  (1), 

Enfin  la  figure  (3)  est  une  representation  dans  le  plan  complexe  de  la  rdponse  d'un  mode  aprAs 
isolation  par  la  mdthode  du  paragraphe  (2).  Elle  reprdsente,  pour  chacun  des  quatre  modes,  le  premier  et 
le  second  membre  de  l'dquation  (2-I7-I)  et  permet  done  d'apprecier  la  qualite  de  l'isolation  modale  donnde 
par  la  mdthode. 

La  comparaison  desbalayages  (figure  2)  et  des  grandeurs  generalises  (tableau  I)  para  ft  satis- 
faisante  si  l'on  dent  compte  du  fait  que  la  structure  essayde  avait  un  comportement  qul  ddnonqait  la 
prdaence  de  non  lindaritds.  La  qualitd  des  rdsultats  est  aensiblement  la  mdme  pour  les  deux  mdthodes,  et, 
dans  ce  cas,  il  semble  que  le  choix  des  utilisateurs  devrait  s'orienter  vers  la  mdthode  de  s.^nthAse  modale 
qui  eat  plus  simple  A  mettre  en  oeuvre  et  plus  rapide  que  la  deuxiAme  mdthode.  Toutefois,  la  separation 
modale  semble  rdserver  des  possibilitds  interessantes  pour  les  structures  comportant  de  faibles  non  linda- 
ritds  :  elle  peut  faciliter  l'analyse  de  ces  non  lindaritds  et  il  existe  un  espoir  de  pouvoir  l'associer  A 
un  mo dA le  de  rdponse  comportant  de  faibles  non  lindaritds. 


CONCLUSION 


Lee  deux  methodes  presentdes  dans  cet  expose  permettent  de  determiner  les  caractdristiques  mo- 
dales  d'une  structure  A  partir  d'un  essai  harmonique  avec  une  excitation  qui  n'a  pas  besoin  d'etre  appro- 
pride  mals  qui  doit,  dvidensaent,  introdulre  suffisamment  d'dnergie  dans  les  modes  que  l'on  veut  identifier. 

A  preniAre  vue,  la  mdthode  de  synthAse  modale  paraft  plus  attrayante  pour  les  applications 
pratiques  parce  qu'elle  est  plus  simple  A  mettre  en  oeuvre  que  la  mdthode  de  separation  modale  et  donne 
des  rdsultats  de  qualitd  sensiblement  dquivalente.  Cependant.  par  son  caractAre  analytique,  qui  rdsulte 
de  l'isolation  des  modes,  la  deuxiAme  mdthode  paraft  plus  apte  A  s'adapter  aux  structures  prdsentant  des 
faibles  non  lindaritda. 
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ANNEXE  A 

PROJECTION  ORTH OGON ALE  D'UNE  FONCTIOM  SUR  UNE  BASE  COMPLEXE 
La  bas*  2  est  la  ligne  del  variable!  d'Atat  ZT  (<u>) 

Soit  une  fonccion  cooplexe  .  On  cherche  le  vecteur  £  qui  rend  minimum  l'Acart  j 

£  =  j  !  Zrc  -  A 

ou  00 

t  =  J(CZ*  A)(ZrC  -  A)  dcu 

-0“ 

si  Q  esc  une  solution,  on  a  :  U  *  o  V  s  c 

est  donnfi  par  :  ^  00 

Se=  Set  z*(,z?c-A)db>  4-  J (cz*-a*)ztj«)  Sc 

~<P  „  -*>  Q0 

S  JC*TJ(Z?2?C  -ztydco  ■*-  Jcy  (zz*TC*-Z  A*)  du> 


La  condition 


-OO  -00 

SB  2.0  V  Sc 


riel  ou  cooplexe  entraine  : 


J(z*ZTC  -  2*a)  dui  2  o 


-so 


d'ofli  : 


00  i  00 

C  a.  [  j  Z*Zr dco]  Jz'aMo 

-00  -OO 


Loraque  la  colonne  ddfinie  par  (1-3)  eat  rfielle,  le  vecteur  d'Atat  Z  donnA  par  (I- 

priAtA  Z  (-*•*)  ,  d'ob  l'on  dAduit  : 

,O0 


La  matrlce 


jWd.  -  J (z“zT+ 2  z’^du 

-<*  n 


Z  oi(«)  eat  done  rAelle  loraque  le  vecteur  excitation  eat  rAel. 


dAflni  par  : 


a-1 


4)  A  la  pro- 


ANNEXE  B 


VECTEURS  PRO  PRES  DU  SYSTEME  FAIBLEtfENT  DISSIPATIF 

On  se  propose  de  determiner  les  vecteurs  vdrifiant  1  'equation 

y[f,-A]  -  c 


*  i-/x  -/*>] 


b-1 


b-2 


ou  M  Ct  Jf  sont  des  matrices  diagonales  et  fi  matrice  symdtrique  dont  les  coefficients  peuvent  fitre 
cons iddrds  comae  infinlment  petits  du  Ier  ordre,  alors  que  les  coefficients  de  JA  et  de  ^  sont  finis. 


4}  ^ 

Considdrons  d'abord  les  solutions  fondamentales  | 6,  r  du  systbme  d'dquations  du  second 

Design  oris  par  Aj  ,  la  diagonale  de  ^3  et  consid£rons  *■>  /3q  conme  une  perturbation  du 
iagonal.  Celui-ci  admet  des  solutions  :  •  " 

(tVtr/WH-0 

Avec  la  partie  non  diagonale  de  Q  les  valeurs  proprea  et  les  vecteurs  propres  deviennent 
respect ivement  :  * 


ordre  : 


systbme  diagonal 


b-5-1 

b-5-2 

b-5-3 

conme  des 


r,  =  r„ 

Nous  admettrons  que  le  vecteur  Ur  est  normd  par  la  condition  : 

urTu  -  i 

qui  est  vdrifide  par  le  vecteur  initial  . 

En  reportant  (b-5-I)  et  (fc-5-2)  dans  (b  -3)  et  en  cons i ddrant  A_rt  >  Ay  et  (1/ 
infiniment  petits  du  Ier  ordre,  on  a  :  I  f®  ”  * r 

Ar  *  T.,(P-fr)er  H)--° 

En  prdmultipliant  (  b-6  )  par  ajf  il  vient  : 

Main  ^3-/3$  ne  contenant  pas  de  terme  diagonal,  0 J  )  ft 

est  nul  et*  que  yj»  s  ^  au  second  ordre  prbs.  '  '  I  *  r 

En  prdmultipliant  (b-6)  par  r  ,  et  en  remarquant  que  fiL*  JA  &  et  /J 

sont  nuls  puisque  et  ^  sont  diagonales,  on  a  :  T  *  I  •> 

Tf  KP  tr*K( 

ou,  puisque  yU^  ^  et  ^  sont  diagonales, 

t,  elpe,  *  < K 

°"r  ,  ,I«H  Pr  Ar  *  4  (  fr'Ak  *RA**ifn)  KW  =  ° 

A  lj/  ,  qui  ddsigne  le  k  dldment  du  vecteur  ^  ,  est  donnd  par  : 

eT|f  * - TrJ*kr - 

On  a  :  As.  s  -  (<fVr  et  l'hypothise  du  syetfcme  faibleoent  dlssipatif  conduit  A  considdrer  cir 

comae  infiniment  petit  du  Ier  ordre,  de  mime  que  Si,  de  plus,  on  atfaet  que  les  pulsations  propres 

sont  eufflswment  sdpardes  pour  que  y>  solt  fini  yjjl  r  vient  : 


b-6 


est  nul.  On  en  ddduit  que 


.  <2 
O  x 
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CT  a/  -  _ 

vr  A(f-ti) 


/ 


oil 


f\  <Wli8Tie  le  tenne  Ak  de  la  matrice  diagonale ( /*ic~  t  ^k  ) 

Le  coefficient  peut  dtre  ddduit  de  la  condition  de  normalisation  (b*5-3)  .  On  a  : 

0,TUr-  teT**S'T)(Vi1'r) 

=  i  (yrer+£'ti) 


d’oO  l’on  ddduit  : 


m  ^  d'apres  (b-5-3) 

En  dSsignant  par  S'  la  matrice  des  vecteura  ^  ,  on  a  done  fina  lenient 

U  =  er  +  -£•  4>er 

r  _  ^  fikf _  ,  W*  r 

er  if/e  = 


b-8 


k  =r 

En  transposant  (b-3)  on  voit,  compte  tenu  de  la  symStrie  des  matrices 


que  Uj  vdrifie  1 'Squat ion  : 


/*,/*  et  r  • 


On  en  dSduit  que  le  vecteur  donnS  par 


Ur  (W’TrP't)'0 

ar  : 

/  m 

*r  =  U/'H 

vdrifie  1’ Squat ion  (b-1)  et  constitue  done  le  vecteur  propre  k  gauche  cherchS  .  On  a  done  : 

*[1*^1  er  \ 


f  er/ 


b-9-l 


oh  eat  la  matrice  ddfinie  par 

<rer  = 


^  (-’kr 


AK-W 


le*r 


War 


b-9-2 


Fig.l  -  Avion  TB  20.  Position  des 
points  de  mesure  et  du 
point  d'excitation 


L 


0 
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*  • 

16 


16 


20 


f)  -  Capteur  733 


l»  ill  an  pr»  maa  | 


Frequence(Hz) 

-  ^  • 

20 


h)  -  Capteur  737 
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Fig .  2  -  Balayage  en  frd' 
quence  antisymdtrique, 
excitation  derive, 
point  728.  Reconstitu¬ 
tion  des  rdponses  Vi¬ 
tesse  . 

-  mesurc 

o  o  synthSse  mo- 
dale  (mdth.  §  1) 

•  •  isolation 
modale  (tndth.  §  2) 


c)-  Mode  3  :  Vibration  masses  d  dqui-  d)-  Mode  4  :  Coplanaire  empennage 

librage  avec  participation  empen-  horizontal 

nage  horizontal  f,»  17,49  Hz,  c<  •  0,041 

ft.”  15,59  Hz,  oC  =  0,012 

Fig. 3  -  Isolation  modale.  Representation  dans  le  plan  complexe  de  la  rdponse  des  modes  isolds 
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Nonlinear  Aeroelastic  Analyses  Taking  Into  Account 
Active  Control  Systems 

by 
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Summary 

Analytical  investigations  on  aircraft  with  active  control  systems  can  be  performed  on  the 
basis  of  an  extended  formulation  of  the  generalized  aeroelastic  equations.  Besides  the 
structural  modal  parameters  needed  for  these  investigations  which  can  be  determined  in  a 
ground  vibration  test  on  the  real  aircraft  structure,  a  comprehensive  knowledge  of  the 
partially  nonlinear  transfer  functions  of  the  different  elements  of  the  control  loops  is 
necessary.  For  this  purpose  a  method  is  presented  which  shows  the  determination  of  the 
transfer  functions  of  these  elements  in  a  ground  vibration  test  on  the  real  aircraft  and 
in  special  tests  on  the  aircraft  hydraulic  actuation  systems.  In  order  to  demonstrate  the 
validity  of  the  elaborated  procedure,  dynamic  response  calculations  based  upon  the  pro¬ 
posed  extended  formulation  of  the  aeroelastic  equations,  which  take  into  account  both  the 
results  of  the  ground  vibration  test  as  well  as  of  the  actuator  tests,  are  performed  on  a 
model  wing  structure  with  a  hydraulically  driven  rudder.  Results  of  these  analytical  in¬ 
vestigations  are  compared  to  measured  data  on  the  same  model  structure. 

1 .  Introduction 

In  recent  years  a  tendency  to  install  increasingly  complex  control  systems  in  modern  air¬ 
craft  can  be  observed,  see  Ref. Ml.  These  control  systems  can  greatly  affect  the  dynamic 
behavior  of  the  aircraft.  Especially  the  integration  of  active  control  systems  operating 
in  the  higher  frequency  range  such  as  active  flutter  suppression  or  gust  alleviation  sys¬ 
tems  (see  Refs. [ 2 , 3 , 4 ] ) ,  which  are  introduced  into  modern  aircraft  design  in  order  to  pro¬ 
voke  a  desired  change  of  the  aircraft's  dynamic  behavior,  requires  detailed  investigations 
for  determination  of  structural  and  control  system  coupling  effects.  Attempting  to  consid¬ 
er  the  dynamic  properties  of  the  active  control  system  in  the  aeroelastic  qualification  of 
the  aircraft  reveals  some  problems  and  requires  solutions  of  the  following  questions: 

a)  What  possibilities  are  given  to  describe  analytically  the  equations  of  motion  of  an 
aircraft  with  active  control  systems  -  which  is  indeed  a  nonconservative  system  - 
in  a  generalized  form,  say  by  means  of  its  eigenmodes  as  elastomechanic  degrees  of 
freedom? 

b)  How  must  the  ground  vibration  test  be  performed  on  the  aircraft  in  order  to  deter¬ 
mine  its  modal  parameters,  such  as  eigenfrequencies,  eigenmodes,  generalized  masses 
and  generalized  damping  coefficients? 

c)  In  what  way  must  the  dynamic  characteristics  of  the  hydraulic  actuators  and  control 
system  feedback  loops  be  determined  and  introduced  into  the  generalized  equations 
of  motion,  such  that  the  dynamic  behavior  of  servo  systems  can  be  fully  considered 
in  aeroelastic  analyses? 

Neither  these  questions  nor  consequential  problems  may  be  regarded  separately  when  seeking 
an  answer.  For  aeroelastic  investigations  it  is  of  great  importance  to  find  one  consis¬ 
tent  approach  which  gives  a  satisfactory  answer  to  all  three  points  mentioned  above. 

In  Ref s . t 5 , 6 , 7 , 8]  several  possible  approaches  are  discussed.  Unfortunately  in  these  reports 
some  of  the  key  problems  such  as  the  determination  of  transfer  functions  of  the  actuators 
or  the  performance  of  ground  vibration  tests  are  not  treated  in  enough  detail  or  in  a 
practical  way.  Furthermore  most  of  these  approaches  are  based  on  the  more  or  less  incorrect 
assumption  of  linear  behavior  of  the  aircraft  with  control  systems. 

This  report  describes  a  consistent  and  practical  new  method  which  can  be  applied  success¬ 
fully  to  the  aeroelastic  qualification  of  aircraft  with  active  control  systems.  First,  in 
Section  2.1,  an  experimental-analytical  method  to  determine  the  frequency  responses  of  the 
actuator  transfer  functions  is  discussed,  taking  into  account  both  nonlinear  and  preload 
effects.  In  Section  2.2  a  description  of  a  ground  vibration  test  procedure  is  given  for 
determination  of  the  modal  aircraft  parameters  as  well  as  structural  and  electric  feedback 
effects  due  to  the  control  system  and  the  servomechanisms.  Finally,  it  is  shown  in  Section 
2.3  how  the  measured  data  resulting  from  the  actuator  tests  (frequency  responses)  and  from 
the  ground  vibration  test  (modal  parameters,  control  system  feedback  properties)  are  intro¬ 
duced  into  the  generalized  equations  of  motion  in  order  to  allow  the  performance  of  aero¬ 
elastic  qualification  analyses  on  aircraft  with  active  control  systems. 

The  aim  of  this  report  is  to  present  a  survey  of  such  procedural  advances.  Further  details 
can  be  taken  from  Ref. (9], 
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2 ,  Description  of  the  method 

Two  of  the  major  problems  encountered  when  performing  aeroelastic  analyses  on  aircraft 
with  active  control  systems  must  be  attributed  to  the  fact  that  such  aircraft  can  be  re¬ 
garded  neither  as  a  linear  nor  as  a  conservative  system.  Due  to  these  nonlinear  and  non¬ 
conservative  properties/  the  limits  of  validity  of  numerous  theoretical  approaches  and 
practical  procedures  commonly  used  for  the  dynamic  qualification  of  aircraft  are  greatly 
exceeded. 

The  nonconservative  as  well  as  significant  (locally  concentrated)  nonlinear  properties 
of  aircraft  with  active  control  systems  must  be  attributed  to  the  actuators.  In  order  to 
allow  the  utilization  of  conventional  ground  vibration  test  procedures,  it  would  be  ad¬ 
vantageous  to  replace  the  (active)  actuators  by  (passive)  dummy  actuators  during  the  per¬ 
formance  of  the  ground  vibration  test.  This  procedure  requires  that  a  possibility  exists 
for  determination  of  the  actuator  transfer  functions  which  must  be  considered  in  the  aero¬ 
elastic  analysis  later  performed  on  the  basis  of  the  modal  parameters  measured  during  the 
ground  vibration  test  on  the  modified  (passive)  aircraft. 

To  satisfy  this  requirement,  a  description  of  a  new  method  for  determination  of  the  non¬ 
linear  and  preload-dependent  frequency  responses  of  the  actuator  transfer  functions  is 
given  in  the  following  section. 

2. 1  An  experimental-analytical  method  for  determination  of  the  frequency  responses  of  the 
transfer  functions  of  hydraulic  actuators 

Numerous  tests  performed  on  actuators  (see  Ref s . [ 1 0 , 1 1 , 1 2 , 1 3) )  have  shown  that  servocylin- 
ders  can  possess  highly  nonlinear,  complex,  frequency-  and  preload-dependent  properties. 

A  further  difficulty  is  caused  by  the  fact  that  actuators  must  be  regarded  as  multivari¬ 
able  systems  with  two  input  signals  and  one  output  signal.  The  input  signals  are  the  servo 
valve  spool  displacement  and  the  actuator  force;  the  output  signal  is  the  actuator  elonga¬ 
tion. 


Fig . 1 :  Rudder/wing  connection  with  hydraulic  actuator. 

The  elaborated  method  for  determination  of  the  frequency  responses  of  the  actuator  trans¬ 
fer  functions  will  be  applied  here  taking  the  example  of  an  actuator  controlled  by  a  me¬ 
chanical  as  well  as  by  an  electric  input  signal  (Fig. 1 ) .  Detailed  analytical  investiga¬ 
tions  described  in  Ref. [9]  have  yielded  that  the  equations  of  motion  of  an  actuator  can 
be  formulated  (in  the  frequency  domain)  either  by  the  force  equation 

(1)  J  =  FE(s)e’  +  FX(s)  (XG-  xK> 

or  by  means  of  the  displacement  equation 

(2)  XQ-  Xy.  =  XE(s)e'  +XF(s)J 

when  making  use  of  the  following  notations  (Gothic  letters  refer  to  Laplace  transformed 
variables) : 

3^.,  XR  :  displacement  of  the  power  cylinder  block  and  of  the  power  piston  rod,  respectively, 
J  :  actuator  force. 
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(>’  :  relative  displacement  of  the  servo  valve  spool  defined  as 


(3) 

e’  «  e  -  xG 

s=o+iw 

:  Laplace  variable, 

FE(s) 

:  transfer  function  relating 
force , 

the 

FX(s) 

:  transfer  function  relating 
(actuator  admittance) , 

the 

XE(s) 

:  transfer  function  relating 
elongation, 

the 

XF  (s) 

:  transfer  function  relating 
(actuator  impedance) . 

the 

servo  valve  spool  displacement  to  the  actuator 
actuator  elongation  to  the  actuator  force 
servo  valve  spool  displacement  to  the  actuator 
actuator  force  to  the  actuator  elongation 


Between  the  actuator  impedance  and  admittance  transfer  functions  exists  the  following  re¬ 
lation: 

(4)  XF (s)  =  IFX(s))'1 

Eqs.(1)  and  (2)  are  valid  if  the  mass  forces  of  some  actuator  components  (piston,  piston 
rod,  power  cylinder  block)  are  neglected.  The  transfer  functions  of  these  equations  are 
assumed  to  have  amplitude-  and  preload-dependent  properties.  It  could  be  shown  that  if  the 
nonlinear  parameters  (flow,  leakage,  damping  coef f icients)  included  in  these  transfer  func¬ 
tions  are  considered  on  the  basis  of  the  "harmonic  balance"  procedure  (see  Ref. [14])  - 
thus  allowing  the  performance  of  a  nonlinear  analysis  -  then  the  following  relationship 
exists  between  the  frequency  responses  of  the  transfer  functions  XE|s),  FE(s)  and  FX(s): 

(5)  XE(iuQ,e^)  =  -  FE(iw0,e^,f0)  •  (FX (i“0 , f g) ) ~ 1 

Performing  this  multiplication  of  frequency  responses  at  discrete  circular  frequencies 
u)Q  requires  an  exact  consideration  of  the  amplitude  levels  of  el  and  fQ  . 

If  the  cylinder  mass  forces  are  introduced  into  the  mathematical  model,  the  forces  acting 
on  the  power  cylinder  block  and  on  the  power  piston  rod,  designated  by  fG  and  fK  ,  re¬ 
spectively,  are  no  longer  equal.  Eqs.(l)  and  (2)  can  then  be  written  as  follows: 

(6a)  =  FKE(s)e'  +  FKXK(S)3CK  -  FkXg(s)Xg  . 

(6b)  fG  =  FqE ( s ) C  +  FgXk1s,3Ck  +  FGXG*S*3CG  ' 

and 

(7)  =  XE<S>  *  +  xgfg<s>  •  fG  +  XGFG,S)  *  mGsJ  ‘  XK 

Denoting  the  mass  of  the  power  cylinder  rod  (with  the  piston)  by  mR  and  the  mass  of  the 
power  cylinder  block  by  itiq  ,  it  can  be  shown  that  with 


(8)  xgfq  <  s )  =  (FgXg<s))-' 

the  following  relations  exist  between  the  transfer  functions  of  Eqs.(6)  and  (7)  and  those 


Of  Eqs 

.  ( 1 )  and 

(2) 

: 

(9a) 

FrE(s) 

= 

FGE<S) 

=  FE (s) 

(9b) 

fkxg(s) 

- 

FX  (s) 

r 

(9c) 

FgXr(s) 

= 

-  FX(S) 

f 

(9d) 

FkXk(s) 

= 

-  FX (s) 

+  mK  *  52 

(9e) 

fgXg(s) 

= 

FX (s)  - 

mG  •  s! 

Considering  the  mass  forces,  Eq. (5)  can  also  be  rewritten  as  follows: 
(10)  XE(i«i0,e’l  =  FGE(iui0,e^,f0)  •  (FQXK  ( ieQ ,  f  Q)  )  _1 


Eqs. (7)  and  (8)  can  be  regarded  as  the  equations  of  motion  of  a  "free-free”  supported 
actuator.  In  the  case  of  the  actuator  drawn  in  Fig.1  the  following  relation  can  be  written 
for  the  relative  displacement  e'  of  the  servo  valve  spool  if  kinematics  of  the  control 
rods  and  a  flexibility  at  the  power  piston  rod  bearing  are  considered: 

(ID  e'  =  -s^s  +  62eSK  +  8  3**-^ 

with 


displacement  of  the  control  stick  rod, 
displacement  of  the  servo  piston, 

lever  arm  rate  of  the  control  rods  as  defined  in  Fig.1. 


Denoting  the  nonlinear  transfer  function  of  the  amplifier  unit  comprising  the  electrohy- 
draulic  servo  valve  and  the  servo  piston  (without  regarding  the  servo  piston  displacement 


feedback  loop) ,  the  servo  piston  displacement  eSK  can  be  related  to  the  electric  input 
signal  6e  of  the  electrohydraulic  servo  valve  as  follows: 

(12>  CSK  *  Gem(s>  «e 


(14)  eeSK  =  Ge(s)  ,eSK  “  XG* 

Thereby  the  following  notations  were  used: 

:  electric  servo  valve  input  signal  induced  by  the  active  control  system, 
eeSK  :  electric  servo  valve  input  signal  induced  by  the  servo  piston  feedback  loop, 
Ge(s)  :  transfer  function  of  the  servo  piston  displacement  feedback  loop. 

Introducing  Eqs. ( 1 2 , 1 3 , 1 4)  into  Eq .(11)  yields: 

(15)  c  =  G«m(s)eI  ♦  eB  *  ef  -  p  ♦  g^og^.))  .  xG 
with 

,16)  Gem(s>  “  *2  ~  G  gTs)G  (s) 

em  e 

as  the  transfer  function  of  the  amplifier  unit  consisting  of  the  electrohydraulic  servo 
valve  and  the  servo  piston  cylinder,  taking  into  account  the  servo  piston  feedback  loop. 
The  servo  valve  spool  inputs  and  due  to  the  displacement  of  the  control  stick 

rod  and  to  the  flexibility  of  the  power  piston  rod  bearing  are  defined  in  accordance  wit) 
Eq.(ll)  as  follows: 

07)  eK  -  -«,xs 


Introducing  Eq.(15)  into  Eqs. (6)  and  (7)  yields  the  equations  of  motion  of  the  most  com¬ 
plicated  actuator  type  controlled  by  a  mechanical  as  well  as  by  an  electric  input  signal. 
Because  the  displacement  equation  of  motion  is  of  primary  importance  for  the  aeroelastic 
analysis  later  performed  in  Section  2.3,  this  equation  is  rewritten  here  in  its  entirety 
yielding: 


XG"XK 


=  xe (s)  {c4(s)  #I  ♦  era  ♦  ef  -  (i  ♦  <4<s)-Ge(s)) -x^}  ♦  ♦  x^G(s). 


"fcs-XK 


Having  obtained  this  equation  from  analytical  investigations,  the  next  problem  to  be  solved 
is  how  to  determine  the  frequency  responses  of  the  transfer  functions  XE(s),  XgFq(s), 

G?m(s)  and  Ge(s).  These  frequency  responses  can  hardly  be  determined  analytically  because 
it  is  very  difficult  to  give  an  accurate  evaluation  of  some  actuator  parameters  such  as 
the  servo  valve  flow,  various  pressure  gradients  or  leakage  coefficients  for  the 
dynamic  operation  of  the  actuator.  It  could  be  shown  (see  Ref.  [15])  that  these  para¬ 
meters  are  extensively  frequency-dependent  and  exhibit  nonlinear  properties.  Thus  the  fre¬ 
quency  responses  of  the  various  transfer  functions  must  be  determined  experimentally  in 
order  to  obtain  reliable  results.  With  the  exception  of  the  frequency  responses  G|m(iu) 
and  Ge(im)  which  can  be  determined  easily  by  a  test,  it  would  be  very  cumbersome  to  meas¬ 
ure  the  (preload-dependent)  frequency  response  XE(iu>)  directly  in  a  test. 

The  key  idea  of  the  method  proposed  here  is  that  frequency  responses  of  transfer  functions 
which  cannot  easily  be  determined  directly  in  a  test  are  calculated  from  frequency  responses 
of  other  transfer  functions  the  experimental  determination  of  which  is  quite  simple.  A  de¬ 
scription  of  this  method  to  determine  the  frequency  responses  of  the  transfer  functions 
listed  in  Eq. (19)  can  be  summarized  as  follows: 

a)  Measurement  of  the  frequency  response  F^X^  ( ius )  ,  (Sec.  2. 1.1), 

b)  Measurement  of  the  frequency  response  FGE(iu>),  (Sec.  2. 1.2), 

c)  Measurement  of  the  frequency  responses  G|m<iui)  and  Ge(iui),  (Sec. 2. 1.3), 

d)  Calculation  of  the  frequency  response  FGXG(iu)  on  the  basis  of  the  measured 

frequency  response  FGXK(iu)  according  to  Eqs. (9c)  and  (9e) . 

Calculation  of  the  frequency  response  XE(iui)  on  the  basis  of  the  measured 
frequency  responses  FGE(iw)  and  FGX^ (iw)  according  to  Eq.(10). 

A  detailed  description  of  the  performance  of  the  various  actuator  tests  will  be  given  in 
the  following  sections.  Because  all  transfer  functions  (except  Ge(s))  are  dependent  on 
the  static  actuator  preload  force,  these  tests  are  performed  at  the  various  levels  of  the 
preload  force  of  interest.  Multiplication  of  frequency  responses  -  as  required  by  Eq .(10)  - 
is  consistent  only  if  all  frequency  responses  are  attributed  to  the  same  preload.  Similar¬ 
ly  two  other  working-point  parameters  have  to  be  considered.  Firstly,  the  frequency  re¬ 
sponses  are  dependent  on  the  static  elongation  of  the  actuator  (static  position  of  the 


power  piston  in  the  cylinder  block) .  Secondly,  since  the  bulk  modulus  of  the  oil  is  in¬ 
fluenced  by  temperature  changes,  frequency  responses  will  also  be  influenced  by  changes 
of  oil  temperature.  In  order  to  avoid  the  influence  of  the  latter  working-point  parame¬ 
ters,  all  following  tests  will  be  performed  at  a  nearly  constant  oil  temperature  and  at 
a  given  static  actuator  elongation. 

2.1.1  Experimental  determination  of  the  frequency  response  of  the  transfer  function 
fgxk (s> 

The  actuator  mounted  in  the  test  rig  shown  in  Figs. 2  and  3  is  excited  at  its  piston  rod 
by  a  power  cylinder.  The  excitation  is  realized  at  various  frequencies  and  amplitudes  of 
the  exciting  force.  During  these  investigations,  the  servo  valve  spool  is  clamped  to  the 
actuator  power  cylinder  block  ( e ’ =  0 )  in  such  a  position  to  maintain  a  desired  static  pre¬ 
load  force  F  .  For  a  given  harmonic  excitation  of  the  actuator  power  piston  rod,  the 
force  fG  (measured  by  means  of  a  load  cell  fixed  at  the  actuator  power  cylinder  block) 
and  the  relative  displacement  xg-xk  between  the  supporting  points  of  the  actuator  (meas¬ 
ured  by  means  of  a  displacement  pickup)  are  determined.  With  the  help  of  co-quad  analysers 
the  first  harmonic  oscillations  are  filtered  out  of  the  measured  periodic  signals.  Simul¬ 
taneously  the  amplitude  and  the  phase  lag  (relative  to  the  harmonic  output  signal  of  the 
frequency  generator)  of  the  respective  first  harmonic  oscillation  signals  are  determined. 


Assuming  that  the  dynamic  displacement  of  the  actuator  power  cylinder  block  is  very  small 
(xq-O)  ,  Eq.  (6b)  yields  for  a  given  circular  frequency  u>q  and  for  a  given  force  amplitude 
fG=f0  with  regard  to  xG=0  and  e'=0: 

(20)  FGXK(i*0,f0)  =  fG/xR 

Therewith  a  point-by-point  determination  of  the  nonlinear  frequency  response  curves 
FGXK(iw)  is  possible.  Consequently  a  plot  of  measured  frequency  response  curves  is  shown 
in  Fig. 4  (dotted  letters  mark  effective  values) .  In  Fig . 5  the  preload  dependency  of  the 
same  actuator  is  plotted  for  a  constant  frequency. 

2.1.2  Experimental  determination  of  the  frequency  response  of  the  transfer  function 
FGE(s)  “  " 

The  actuator  is  mounted  "rigidly"  in  the  test  rig  shown  in  Figs .6  and  7.  During  the  per¬ 
formance  of  the  test,  the  servo  piston  is  clamped  to  the  actuator  power  cylinder  block 
(esK~*G=0) .  Using  an  (electrodynamic)  exciter,  the  control  rod  is  harmonically  excited  at 
various  displacement  amplitudes  and  frequencies.  For  a  given  excitation  the  relative  dis¬ 
placement  e'  of  the  servo  valve  spool  (by  means  of  a  displacement  pickup)  as  well  as 
the  actuator  force  fG  (by  means  of  a  load  cell)  are  measured. 


Neglecting  the  flexibilities  of  the  test  rig  (xG=xK=0)  ,  Eq. (6b)  yields  for  a  given  cir 
cular  frequency  wq  and  for  a  given  exciting  amplitude  e'se^el  with  respect  to 
xG=xK=0  and  ej=0: 


<21»  rGE<1“0'a0»  *  Ve< 

Consequently  Fig. 8  shows  the  frequency  response  plot  of  the  transfer  function  FQE(s) 
measured  on  an  F-4F  high  gain  aileron  actuator. 


Ihspioctmeni 
Pickup  ■> 
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Fig. 8;  Frequency  response  plot  of  the 
transfer  function  FqE(s) 

(F-4F  aileron  actuator) . 


Fig. 9:  Test  setup  for  determination  of  the 
frequency  responses  GR (lw)  and 
Ge(iu>). 


2.1.3  Experimental  determination  of  the  freguencv  responses  of  the  transfer  functions 


GR  (s)  and  G_(s) 
em  e 


While  performing  this  test,  the  actuator  remains  in  the  test  rig  shown  in  Figs. 6  and  7. 
During  the  test  the  control  rod  is  clamped  to  the  actuator  power  cylinder  block  (em=0) . 

The  electrohydraulic  servo  valve  is  then  harmonically  excited  at  various  frequencies  and 
at  various  amplitudes  of  the  input  voltage  ej  (see  Fig ■ 9 ) .  For  a  given  excitation  ej  , 
the  relative  displacement  e'  of  the  servo  valve  spool,  the  relative  displacement 
(eSv-xG)  of  the  servo  piston  as  well  as  the  electric  signal  eegn  induced  by  the  servo 
piston  displacement  feedback  loop  are  measured.  For  a  given  circular  frequency  u>0  and  a 
given  amplitude  of  the  input  voltage  ei=ejQ  ,  Eq. (15)  yields  with  e^O  ,  ef=0  and 
xG=°: 

(22)  G«m(iUo,eI0)  =  e’/ej 

Furthermore  with  the  assumption  of  a  linear  transfer  behavior  of  the  servo  piston  displace¬ 
ment  feedback  loop,  Eq.(14)  yields: 


Ge(iu>0) 


“eSK^ (eSK 


Having  terminated  the  tests  described  in  Sections  2.1.1,  2.1.2  and  2.1.3,  all  data  are 
available  to  calculate  the  frequency  responses  of  the  transfer  functions  of  the  actuator 
displacement  equation  of  motion  (£q.(19)). 

It  must  be  emphasized  that  all  these  investigations  on  actuators  were  carried  out  with  the 
Intention  to  replace  the  active  actuators  by  passive  dummy  actuators  during  the  performance 
of  the  ground  vibration  test  on  the  real  aircraft  structure.  The  requirement  formulated  at 
the  beginning  of  Chapter  2,  which  stated  that  the  replacement  of  the  actuators  by  dummy 
actuators  -  leading  to  a  modification  of  the  aircraft  -  is  practicable  only  if  determina¬ 
tion  of  the  actuator  transfer  functions  is  possible,  is  now  satisfied.  Finally  it  should 
be  mentioned  that  actuator  tests  as  described  in  Sections  2.1.1  through  2.1.3  could  hardly 
be  performed  with  such  accuracy  and  at  such  high  preload  and  dynamic  exciting  force  levels 
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in  a  ground  vibration  test  on  the  aircraft  equipped  with  the  actuators.  Having  solved  the 
actuator  problem,  attention  can  now  be  focussed  on  the  ground  vibration  test  procedure. 

2.2  Performance  of  the  ground  vibration  test  on  aircraft  with  active  control  systems 

For  the  performance  of  aeroelastic  analyses  on  the  basis  of  modal  data  it  is  necessary 
that  the  eigenmodes  fulfill  the  orthogonality  condition.  Practically  speaking,  determina¬ 
tion  of  orthogonal  eigenmodes  in  the  ground  vibration  test  is  possible  only  if  the  struc¬ 
ture  has  conservative  and  nearly  linear  properties. 

Dealing  with  structural  nonlinearities,  a  distinction  must  be  made  between  the  continuous¬ 
ly  distributed  (global)  structural  nonlinearities,  the  origin  of  which  is  located  in  the 
riveted,  screwed  or  bolted  connections,  and  the  locally  concentrated  nonlinearities  due 
to  components  with  nonlinear  properties  such  as  actuators  or  partially  filled  external 
store  tanks.  In  Ref. [16]  it  is  shown  that  weak  global  nonlinearities  are  of  primary  im¬ 
portance  only  when  dealing  with  neighboring  eigenf requencies .  Locally  concentrated  non- 
linearities  located  in  the  rudder/wing  connections  or  in  external  store  systems  can  be 
eliminated  in  the  ground  vibration  test  by  replacing  the  nonlinear  elements  (actuators, 
stores)  by  dummy  components.  A  subsequent  consideration  of  the  true  nonlinear  behavior  of 
the  replaced  components  in  the  generalized  equations  of  motion  is  possible  by  adding  mod¬ 
ification  terms  in  these  equations  (see  Refs. [17, 18, 19] ) . 

The  consequence  of  replacing  the  active  actuators  with  passive  dummy  actuators  having  lin¬ 
ear  stiffness  characteristics  is  that  such  modified  aircraft  can  be  regarded  as  an  elasto- 
mechanical  structure  with  quasi-linear  and  conservative  properties.  This  fact  allows  appli¬ 
cation  of  conventional  methods  for  the  performance  of  the  ground  vibration  test. 

2.2.1  Determination  of  the  dynamic  properties  of  the  rudder  rigid  body  degrees  of  freedom 

Since  the  aeroelastic  analyses  are  performed  later  in  Section  2.3  on  the  basis  of  quasi¬ 
normal  coordinates  (see  Ref. [20]),  characterized  by  the  fact  that  rudder  rigid  body  modes 
and  global  aircraft  eigenmodes  are  considered  separately  in  the  generalized  equations  of 
motion,  it  is  advantageous  to  determine  the  modal  parameters  of  global  aircraft  eigenmodes 
and  rudder  rigid  body  modes  in  separate  tests  as  well. 


Fig. 10:  Rudder/wing  connection  with  dummy  actuator. 


Before  the  performance  of  the  ground  vibration  test,  the  hydraulic  actuators  are  replaced 
by  dummy  actuators  (Fig .10).  Since  Eq.(l9)  requires  the  mechanical  servo  input  signals 
em  and  e^  and  the  electric  signal  ej  to  be  determined,  a  displacement  pickup  (located 
in  the  dummy  servo  valve)  and  a  load  cell  measuring  the  actuator  force  are  fitted  in  the 
dummy  actuator.  The  measurement  of  the  electric  input  signal  ej  can  be  achieved  with  the 
help  of  a  voltameter.  Furthermore  for  determination  of  the  actuator  attachment  stiffnesses, 
an  acceleration  pickup  is  placed  on  the  dummy  actuator  rudder  bearing. 

Before  mounting  the  various  dummy  actuators  in  the  respective  rudder/wing  connections, 
their  stiffness  characteristics  are  determined  in  preliminary  tension  tests.  The  stiff¬ 
nesses  of  the  dummy  actuators  are  chosen  such  that  the  eigenfrequencies  of  the  rudder 
rigid  body  modes  do  not  fall  into  the  (immediate)  frequency  range  of  a  global  aircraft 
eigenmode.  To  satisfy  this  condition,  resultB  of  previously  performed  (finite  element) 
calculations  have  to  be  considered. 

The  connecting  parts  of  the  dummy  actuators  have  to  be  constructed  within  the  same  manu¬ 
facturing  tolerances  as  for  the  original  actuators.  If  the  original  connecting  parts  are 
fixed  by  screw  connections  on  the  original  actuators,  it  will  be  possible  in  most  cases 
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to  use  these  original  fittings  for  the  fixation  of  the  dummy  actuators.  Also  the  original 
connection  parts  such  as  bolts  or  lugs  will  be  used  when  mounting  the  dummy  actuators  in 
the  rudder/wing  connections. 

Moreover  the  freedom  of  movement  of  the  steering  mechanisms  {control  stick,  foot  pedals) 
is  reduced  by  attaching  these  mechanisms  to  the  airframe  with  (rubber)  cords.  Thus  unde¬ 
sirable  movement  of  the  control  stick  and  foot  pedal  due  to  friction  effects  in  the  bear¬ 
ings  of  the  steering  mechanism  rods  can  be  avoided  when  the  aircraft  is  excited. 

While  performing  the  ground  vibration  test,  the  dynamic  properties  of  the  rudder  rigid 
body  degrees  of  freedom  are  determined  one  after  another.  Thereby  the  p-th  rudder 
<p=1,2,...m)  is  excited  in  its  eigenfrequency  by  a  (generalized)  force  Rp  .  Considering 
the  notation  used  in  the  mathematical  model  of  Fig. 10  as  well  as  the  block  diagram  of 
Fig. 1 1  yields: 

>24>  Rp  "  "p  '  Bp 

with 


R^  :  generalized  external  rudder  exciting  force, 
n^  :  rudder  hinge  moment, 

Bm  :  rudder  hinge  angle. 


Fig. 1 1 :  Block  diagram  showing  the  coupling  between  the  r-th  global  aircraft 
eigenmode  and  the  u-th  rudder  degree  of  freedom  of  the  modified 
aircraft  with  dummy  actuators  (ground  vibration  test  configuration) . 


Assuming  that  none  of  the  global  aircraft  eigenmodes  are  excited  by  this  force  -  since 
there  is  no  immediate  eigenfrequency  neighborhood  of  rudder  modes  and  global  aircraft  ei¬ 
genmodes  -  the  p-th  rudder  rigid  body  mode  can  be  regarded  as  a  discrete  degree  of  free- 


dom, 
(r=1 , 

the 

2,.. 

equations  of  motion  of  which  can  be  formulated  as  follows  with 
.n)  : 

(25) 

fu 

= 

■  C1U  XKu 

(26) 

£u 

X 

CDtl(XKp  "  XGu’  ' 

(27) 

nu 

= 

‘  fu  YU 

(28) 

with 

* 

XGu/yu  +  n(/C2p 

qr  :  generalized  coordinate  of  the  r-th  global  aircraft  elgenmode, 
f„  :  dummy  actuator  force. 

Cup  s  dummy  actuator  stiffness, 

Yw  :  rudder  lever  arm, 

C1(J  :  nonlinear  global  attachment  stiffness  on  the  wing  side, 

C2W  :  nonlinear  global  attachment  stiffness  (in  torsion)  on  the  rudder  side. 
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Solving  'Eqs.  (25)  to  (28)  with  respect  to  and  C2p 


(29)  Clu(fu) 


fU  *  XGu  CDy 


(30)  C,„ (f„) 


IbJjl. 


24  u  ■  x  /y  =  B 

Gu  4  4 


yields : 


Utilizing  the  pickups  mounted  on  the  dummy  actuator,  the  actuator  force  and  the  dis¬ 

placement  Xqj,  are  measured  at  various  exciting  force  levels.  The  hinge  angle  Bu  is 
also  measured.  Introducing  the  measured  data  into  Eqs.  (29)  and  (30)  allows  the  nonlinear 
attachment  stiffness  characteristics  to  be  determined.  If  the  stiffness  parameters  are 
available  from  static  tests,  these  investigations  are  optional  or  can  be  used  in  order  to 
check  static  test  data. 


Performing  the  test  mentioned  above,  the  servo  valve  spool  displacement  ep=efp  -  due  to 
the  flexibility  of  the  actuator  attachment  on  the  wing  side  -  is  also  measured  by  means  of 
the  servo  valve  displacement  pickup.  The  nonlinear  transfer  function  of  structural  recou¬ 
pling  hC(j  can  then  be  calculated  as  follows: 


(31) 


heu(V 


=  Cf4/f4 


Further  investigations  are  performed  to  determine  the  generalized  masses  and  generalized 
damping  coefficients  of  the  rudder  rigid  body  modes.  Thus  after  termination  of  these  tests 
the  following  data  are  available: 


•  the  rudder  rigid  body  modes  with  respect  to  the  corresponding  normalized  unit 

hinge  angles  B0m  (h  =  1 , 2 , . . . ,m) , 

•  the  generalized  masses  , 

•  the  generalized  damping  coefficients  , 

•  the  nonlinear  attachment  stiffnesses  C1(J(fp)  and  C2M(fp)  ' 

•  the  nonlinear  transfer  function  of  structural  recoupling  h^lf^)  . 


All  parameters  made  available  by  these  tests  can  be  introduced  directly  into  the  aeroelas- 
tic  analysis  performed  in  Section  2.3. 

2.2.2  Determination  of  the  modal  parameters  of  the  global  aircraft  eigenmodes 

During  the  performance  of  the  tests  described  in  this  section,  the  electric  part  of  the 
active  control  system  is  in  operation.  The  determination  of  the  modal  parameters  (eigen- 
frequencies  v r  ,  orthoqonal  global  aircraft  eigenmodes  $£  ,  generalized  masses  Mrr 
and  generalized  damping  coefficients  Drr)  is  achieved  in  a  conventional  ground  vibration 
test  (classical  phase  resonance  method,  see  Ref. [21]).  It  will  be  assumed  here  that  due 
to  the  eigenfrequency  spacing  between  rudder  rigid  body  modes  and  global  aircraft  eigen¬ 
modes,  there  are  no  rudder  rigid  body  motions  contained  in  the  various  global  aircraft 
eigenmodes.  If  this  assumption  is  not  fulfilled,  the  measured  ground  vibration  test  data 
must  be  corrected  by  a  modification  analyses  (see  Ref s . [9 , 22 , 23] ) . 


Besides  the  modal  parameters  mentioned  above,  the  mechanical  servo  valve  spool  displace¬ 
ments  eur  and  the  electric  servo  valve  (feedback)  input  signals  Ceur  are  determined 
in  each  eigenmode.  The  double  indexing  of  terms,  e.g.  the  force  fyr  ,  means  that  this 
term  is  related  to  the  4-th  actuator  and  to  the  r-th  eigenmode.  Considering  the  block  dia¬ 
gram  of  Fig.  11,  it  can  be  seen  that  the  mechanical  input  signal  e^j-  is  composed  as  fol¬ 
lows  out  of  2  different  parts: 


(32) 


ur 


e  _  +  e  , 

04r  fur 


with 


e0ur  :  mechanical  structural  feedback  signal  at  zero  actuator  force, 

efur  !  mechanlcal  structural  feedback  signal  due  to  the  flexibility  of  the  actuator 
attachment  on  the  wing  side. 


According  to  Eq.(31)  the  following  relation  can  be  written  for  the  second  term  in  Eq.(32): 


(33) 


=  h_ 


"fur  "eu  ur 
Introducing  Eq. (33)  into  Eq. (32)  yields: 


(34) 


e0ur 


ur 


-  h 


eu 


ur 


With  the  assumption  of  linear  behavior  of  the  transfer  functions  characterizing  structural 
feedback  effects  via  the  steering  control  systems,  the  coefficients  Hur  of  the  mechanical 
structural  feedback  matrix  M  are  defined  as  follows: 


(35) 


ur 


e0ur 


xQ/x 


rmax 


with  X0/xrni£x  as  a  normalization  factor.  This  factor  has  to  be  attributed  to  the  normal¬ 
ization  of  tne  various  eigenmodes  by  assigning  the  value  x«  (generally  xg*1mm)  to  the 
measured  amplitude  xrmax  of  the  reference  point  of  the  r-th  eigenmode.  The  (complex) 
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coefficients  Hjjr  of  the  matrix  He  of  electric  (active  control  system)  feedback  are  de¬ 
fined  in  the  same  way  by  the  relation 

(36)  H®r  =  ee(ir  *  xo,/xrmax 

It  has  to  be  postulated  that  due  to  the  filtering  and/or  phase  shifting  properties  of  the 
electric  feedback  systems,  the  coefficients  H®r  can  be  extensively  frequency-dependent. 
The  frequency-dependent  behavior  of  the  electric  part  of  the  active  control  system  cannot 
be  determined  in  the  ground  vibration  test,  for  an  excitation  of  the  structure  in  its  ei- 
genmodes  is  only  possible  at  those  frequencies  corresponding  to  the  respective  eigenfre- 
quencies.  The  frequency  dependence  of  the  coefficients  H[jr  has  to  be  determined  either 
by  calculation  or  in  additional  tests  as  will  be  shown  in  Sec. 2. 2. 3. 

After  completion  of  the  tests  performed  in  this  section,  the  following  data  are  available: 

•  the  eigenfrequencies  vr  of  the  global  aircraft  eigenmodes, 

•  the  global  aircraft  eigenmodes  , 

•  the  generalized  masses  Mrr  , 

•  the  generalized  damping  coefficients  Drr  , 

•  the  coefficients  H,jr  of  the  matrix  of  mechanical  structural  feedback, 

•  the  coefficients  H®_  of  the  matrix  of  electric  feedback. 


2.2.3  Determination  of  frequency  responses  of  the  matrix  of  electric  feedback  coefficients 

As  mentioned  in  Sec. 2. 2. 2  the  determination  of  the  frequency-dependent  properties  of  the 
matrix  of  electric  feedback  in  the  ground  vibration  test  is  impossible.  The  frequency  de¬ 
pendence  of  the  electric  feedback  systems  can  be  determined  by  tests  performed  on  the  ac¬ 
tive  control  system  mounted  outside  the  aircraft.  The  best  testing  installation  for  these 
investigations  is  the  "iron  bird"  test  setup. 


On  this  setup  calibration  tests  are  run  for  each  sensor  by  exciting  the  sensors  one  after 
another  at  a  given  amplitude  y^  and  at  different  frequencies.  Exciting  the  C-th  sensor 
(£  =  1, 2,... P)  the  input  voltage  eefl^(iui)  is  measured  at  all  n  electrohydraulic  servo 
valves  (u=1 ,2, . . .m) .  Knowing  the  normalized  amplitude  of  oscillation  x^r  of  the  i-th 
sensor  in  the  r-th  eigenmode,  the  frequency-dependent  transfer  function  H®r(i&;)  can  be 
calculated  as  follows  when  assuming  linear  behavior  of  the  electric  part  of  the  active 
control  system: 

P 

(37)  H®r(iu.)  =  2^  ee(l5(iu)  •  x5r/y^ 

In  the  case  that  all  transfer  functions  of  the  various  links  of  the  electric  feedback 
loops  (compensation  network)  as  well  as  of  the  sensors  are  known,  determination  of  the 
frequency  responses  Hgr(iu)  is  possible  by  pure  calculation. 

A  check  of  the  frequency  responses  H§r(iaj)  determined  here  is  possible  by  comparing  these 
responses,  formed  at  the  frequencies  corresponding  to  the  aircraft  eigenfrequencies  of  the 
global  aircraft  eigenmodes,  to  the  corresponding  coefficients  Hjjr  measured  in  the  ground 
vibration  test.  By  the  end  of  this  chapter  all  parameters  needed  for  the  performance  of 
aeroelastic  analyses  on  aircraft  equipped  with  active  control  systems  have  been  determined. 


2. 3  Consideration  of  active  control  systems  in  the  generalized  aeroelastic  equations 


Formulation  of  the  aeroelastic  equations  on  the  basis  of  quasi-normal  coordinates  expres¬ 
sing  the  elastic  displacements  of  the  aircraft  by  the  relation 

(38)  X(t)  =  *F  •  q(t)  +  •  p(t) 

with 


*(t)  :  vector  of  the  elastic  displacements  of  the  structure, 

♦  F 

:  modal  matrix  of  the  global  aircraft  eigenmodes, 

♦  R  :  modal  matrix  of  the  rudder  rigid  body  degrees  of  freedom, 
q(t) ,  p(t)  :  vectors  of  the  generalized  (normal)  coordinates, 
leads  to  the  following  generalized  equations  of  motion: 


(39) 


'  Ars  ;  v ' 

/  \ 
{Rr> 

.  v  !  v  _ 

< 

.V. 

>  -  < 

(R  ) 

l  »  ) 

VK  1  0 

f  \ 

{q  ) 

^  1 

r 

(  V 

< 

0  K  „ 

(p  ) 

1  U4  n 

l  “  J 

.A 


i 


30 


with 

M  :  coupled  generalized  mass  matrix, 

D  :  generalized  matrix  of  structural  damping, 

K  :  generalized  stiffness  matrix, 

A  :  matrix  of  the  generalized  unsteady  (motion-induced)  aerodynamic  forces, 
R  :  vector  of  the  generalized  external  forces. 


Assuming  that  the  coefficients  of  the  matrix  A  of  the  generalized  unsteady  aerodynamic 
forces  can  be  determined  either  experimentally  or  theoretically  (see  Ref. [24]),  all  matrix 
coefficients  on  the  left-hand  side  of  Eq. (39)  are  available  from  the  tests  performed  in 
Sections  2.2.1  through  2.2.3  with  exception  of  the  terms  Mru,  M^,  Krr,  K^. 

The  coupling  terms  of  the  symmetric  generalized  mass  matrix  can  be  determined  on  the  basis 
of  a  discrete  mass  matrix  in  of  the  different  rudders,  the  rudder  rigid  body  modes  and 
the  global  aircraft  modes  as  follows: 

(40)  Mr|1  =  Mur  =  (fF)Tm^R  ,  (r=1 , 2 . n  ;  u  =  1,2 . m) 

Furthermore  the  coefficients  Kr-  of  the  generalized  stiffness  matrix  are  related  to  the 
coefficients  Mrr  of  the  generalized  mass  matrix  as  follows  (see  Ref. [9]): 

(41)  Krr  =  wj.  •  Mrr 

with  u>r  as  the  circular  eigenfrequency  of  the  r-th  global  eigenmode. 


Now  the  missing  stiffness  coefficients  of  the  rudder  equations  of  motion  can  be  de¬ 

termined  with  the  help  of  the  actuator  displacement  equation  of  motion  formulated  in  Eq.(l9). 
Taking  this  equation  into  consideration  for  the  u-th  actuator,  we  put: 

<42>  emu  =  *  * 

(43)  el4  =  M^Uw)  •  4 

with  and  MS  as  the  u-th  line  vectors  of  the  matrices  H  and  He  ,  respectively, 

and 

(44)  ef[l  =  hekl  •  fK(i 

The  expression  for  the  total  relative  displacement  of  the  servo  valve  spool  can  then  be 
formulated  as  follows: 


(45) 
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Furthermore  the  following  compatibility  and  equilibrium  conditions  can  be  formulated  (in 
the  same  way  as  in  Section  2.2.1)  for  the  u-th  rudder  rigid  body  degree  of  freedom: 
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(47)  nu 

=  "  fGu  ’  Y“ 

<48>  *Gu 

=  V‘u  -  VC2u> 

Introducing  Eqs. (45)  through  (48)  into  Eq. 
for  the  hinge  angle  of  the  u-th  rudder: 

(19)  finally  yields  the  following  expression 
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If  the  actuator  mass  forces  are  neglected  in  these  equations,  say  mG^mv-O  ,  Eqs.  (49)  and 
(50)  can  be  formulated  as  follows: 
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Although  Eqs. (49)  and  (50)  can  easily  be  handled  in  a  computer  software  program,  further 
considerations  will  be  based  here  on  Eqs. (51)  through  (53).  Introducing  the  following  re¬ 
lations 


with  as  the  normalized  unit  hinge  angle  of  the  u-th  rudder  degree  of  freedom  into 

Eqs. (51)  to  (53)  yields: 
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agonal  matrix  only  if  all  (mechanical  and  electric)  feedback  effects  are  neglected.  Nor¬ 
mally  this  matrix  will  be  filled  in  its  lower  left-hand  part  by  the  coupling  terms  Kur  . 
In  the  case  of  a  harmonic  motion,  Eq.  (39)  can  then  be  rewritten  as  follows: 
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Eq. (59)  describes  the  extended  generalized  aeroelastic  equations  of  aircraft  with  active 
control  systems.  It  has  to  be  noted  that  due  to  the  actuator  and  electric  feedback  loop 
characteristics,  the  coefficients  and  Kyr  are  complex  as  well  as  frequency-  and 

amplitude-dependent.  Calculating  the  dynamic  response  of  the  aircraft  structure  to  a  har¬ 
monic  external  excitation  {Rr,R(i)'r  (the  roofed  letters  mark  peak  pulse  amplitudes),  Eq.  (59) 
has  to  be  solved  with  respect  to  the  vector  of  generalized  coordinates  (qr<Pu)T  ■  Intro¬ 
ducing  this  vector  into  Eq. (38)  permits  the  determination  of  the  elastic  displacements  of 
the  aircraft  structure.  All  mathematical  relations  formulated  in  this  section  are  clearly 
displayed  in  the  block  diagram  of  Fig. 1 2  (the  aerodynamic  forces  are  neglected) . 


Comment:  Since  the  transfer  properties  of  the  actuators  were  determined  in  the  form  of 
frequency  responses,  valid  only  in  the  case  of  a  harmonic  actuator  excitation, 
the  extended  generalized  aeroelastic  equations  were  formulated  here  for  the 
special  case  of  harmonic  aircraft  motion.  Due  to  the  nonlinear  effects  it  is 
not  possible  to  calculate  the  (total)  system  response  to  a  (nonharmonic) 
periodic  external  excitation  by  superposition  of  harmonic  responses  to  different 
external  excitations  (determined  by  a  Fourier-analysls  of  the  periodic  exciting 
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forces) .  A  mathematical  description  of  the  actuator  transfer  functions  in  the 
frequency  domain  (formulation  in  form  of  the  Laplace  variable  s  )  is  possible 
when  approximating  "representative"  frequency  response  curves  -  which  can  be  re¬ 
garded  as  curves  of  Bode  plots  -  by  asymptotic  curves  yielding  the  different 
time  constants  of  the  transfer  functions  (see  Ref. [25]).  By  this  step  the  per¬ 
formance  of  a  nonharmonic  linear  analysis  in  the  frequency  domain 
would  become  possible. 


13i  Rudder /wing  model  structure. 


Hydraulic  actuator. 


3.  Application  of  the  method  to  the  example  of  a  rudder/wlng  model  structure 


In  order  to  check  the  practicability  and  accuracy  of  the  method  described  in  Chapter  2, 
the  method  was  applied  to  the  example  of  a  wing  model  structure  with  a  rudder  driven  by  a 
hydraulic  actuator  (Figs. 13  and  14) .  The  actuator  used  for  driving  the  rudder  was  equipped 
with  only  an  electrohydraulic  servo  valve,  thus  only  electric  (control)  signals  could  be 
used  as  servo  valve  input  signals. 


Fig. 15:  Frequency  response  plot  of  Fig. 1 6 :  Quasi-static  characteristic  of 

the  transfer  function  F_Ee(s).  the  actuator  force  versus  the 

servo  valve  input  voltage. 


3 . 1  Actuator  tests 

The  actuator  tests  were  performed  as  described  in  Section  2.1.  As  a  result  Fig . 1 5  shows  a 
frequency  response  plot  of  the  measured  transfer  function  FgEe(s)  (comparable  to  the 
transfer  function  FqE(s)  listed  in  Eq. (6b)  in  the  case  of  an  actuator  with  mechanical 
and  electric  controls) .  Due  to  the  nonlinearity  in  the  neighborhood  of  zero  force  of  the 
quasi-statically  determined  force-electric  input  voltage  characteristic  (Fig, 16) ,  the  mag¬ 
nitude  characteristics  of  the  plot  shown  in  Fig. 15  are  in  good  agreement  with  the  linear¬ 
ized  values  taken  from  the  quasi-static  curve  only  at  high  levels  of  the  electric  excita¬ 
tion  (and  at  low  frequencies).  For  stability  reasons  actuators  controlled  by  electric  in¬ 
put  signals  have  to  be  equipped  with  a  displacement  pickup  operating  in  a  feedback  loop 
(Fig. 17a).  In  the  case  of  the  actuator  used  here,  this  displacement  pickup  was  integrated 
in  the  actuator  power  cylinder  block.  During  the  performance  of  the  actuator  tests,  also 
the  frequency  response  of  the  transfer  function  Gx(s)  of  this  displacement  feedback  loop 
was  determined. 

3.2  Ground  vibration  test 

The  ground  vibration  test  was  performed  in  the  frequency  range  from  0  to  140  Hz  where 
6  elgenfrequencies  were  determined.  The  eigenmodes  and  the  modal  parameters  corresponding 
to  the  two  lowest  elgenfrequencies  are  drawn  in  Figs. 1 7b  and  1 7c .  Fig. 17d  shows  the  rudder 
rigid  body  mode.  While  performing  the  tests,  no  attention  had  to  be  payed  to  the  attach¬ 
ment  stiffnesses  Ci  and  C2  »  for  the  values  of  these  stiffnesses  were  by  far  higher 
then  the  actuator  stiffness  taken  into  consideration  in  the  later  analysis. 


b)  First  bending  eigenmode. 

c)  First  torsional  eigenmode. 

d)  Rudder  degree  of  freedom. 


3 . 3  Comparison  between  measured  and  calculated  results 

To  check  the  method  proposed  here,  dynamic  response  tests  and  calculations  were  performed 
on  the  rudder/wing  model  structure.  Thereby  the  wing  was  excited  by  a  swept  sine  signal 
with  the  help  of  an  electrodynamic  exciter  (see  Fig.l7a).  The  aim  of  the  investigations 
was  to  determine  the  rudder  hinge  angle  B  (in  form  of  a  polar  plot)  relative  to  the  ex¬ 
citing  force  signal.  The  dynamic  response  calculations  were  performed  in  accordance  with 
the  procedure  described  in  Section  2.3  on  the  basis  of  the  data  resulting  from  the  actua¬ 
tor  tests  as  well  as  from  the  ground  vibration  test  delineated  in  Sections  3.1  and  3.2, 
respectively. 


Fig. 18:  Polar  plot  of  the  rudder  hinge  Fig. 19:  Polar  plot  of  the  rudder  hinge 

angle  8  (Exciting  force  amplitude  =  8N) .  angle  B  (Exciting  force  amplitude  =  ION). 
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As  a  result  of  these  Investigations  2  polar  plots  of  the  rudder  hinge  angle  B  (related 
to  2  different  levels  of  the  exciting  force  amplitude)  are  drawn  in  Figs. 18  and  19.  Be¬ 
cause  of  stability  reasons  the  displacement  feedback  loop  shown  in  Fig .17a  had  to  be  con¬ 
sidered  when  performing  these  investigations.  Figs. 18  and  19  show  good  agreement  between 


Fig, 20:  Frequency  response  of  the 

acceleration  feedback  loop. 


Fig. 21 :  Polar  plot  of  the  rudder  hinge 
angle  6  (Consideration  of  the 
acceleration  feedback  loop) . 
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Fig. 22:  Polar  plot  of  the  rudder  hinge 
angle  8  (Consideration  of  dif¬ 
ferent  acceleration  feedback  loops) . 


measured  and  calculated  results.  In 
these  figures  the  first  polar  plot 
"circle"  is  attributed  to  the  wing 
first  bending  mode,  the  second  "cir¬ 
cle"  to  the  rudder  degree  of  freedom 
and  the  third  to  the  wing  second 
bending  mode.  The  45-degree  rota¬ 
tion  of  the  first  "circle"  (com¬ 
pared  to  the  case  of  a  real  rudder/ 
wing  attachment  stiffness)  is  due 
to  the  phase  lag  characteristics  of 
the  actuator  admittance. 

In  a  second  step  the  acceleration 
feedback  loop  shown  in  Fig. 17a  was 
added.  For  stability  reasons  a  fil¬ 
tering  element  was  integrated  into 
the  electric  feedback  network.  The 
resulting  frequency  response  of  the 
acceleration  feedback  loop  is  drawn 
in  Fig .  20  .  The  corresponding  polar 
plot  of  the  rudder  hinge  angle  B 
is  shown  in  Fig . 2 1 .  Good  agreement 
can  be  observed  between  measured 
and  calculated  results  in  the  case 
of  the  wing  first  bending  mode  and 
the  rudder  rigid  body  mode.  The 
discrepancy  in  the  wing  second  bend¬ 
ing  mode  case  is  due  to  the  effect 
of  actuator  saturation.  Actuator 
data  valid  for  the  saturation  con¬ 
dition  could  not  be  introduced  into 
the  calculations,  for  the  previous 
actuator  tests  had  not  been  per¬ 
formed  at  such  high  actuator  force 
levels.  Fig. 21  shows  that  due  to 
the  acceleration  feedback  loop, 
there  is  a  further  rotation  of  the 
first  polar  plot  "circle"  of  about 
20  degrees.  The  phase  lag  increase 
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has  to  be  attributed  to  the  filtering  element  in  the  acceleration  feedback  loop  and  to 
the  phase  lag  characteristics  of  the  frequency  response  XEe(i.»)  .  The  rather  oval  form 
of  the  polar  plot  “circles"  is  due  to  the  nonlinear  properties  of  the  frequency  response 
XEeiiw). 

The  good  agreement  between  measured  and  calculated  results  shown  in  Figs. 18, 19  and  21  con¬ 
firms  the  consistency  of  the  method  elaborated  here.  Finally  an  interesting  result  may  be 
shown,  obtained  by  pure  calculation  on  the  rudder/wing  model  structure.  Fig. 22  delineates 
the  effect  of  a  phase  change  in  the  acceleration  feedback  loop  on  the  rudder  hinge  angle 
polar  plot.  Therefrom  it  is  obviously  of  paramount  importance  to  consider  the  electric 
(active  control  system)  and  mechanical  (deformations  of  the  steering  mechanism  rods)  prop¬ 
erties  of  the  feedback  systems  in  a  correct  way  when  performing  aeroelastic  qualification 
analyses  and  to  exercise  care  when  determining  the  dynamic  properties  of  these  systems, 
for  they  can  largely  affect  the  dynamic  behavior  of  the  aircraft. 

4 .  Conclusion 

Introducing  active  control  systems  into  modern  aircraft  design  requires  the  development 
of  special  test  and  analytical  procedures  for  the  aeroelastic  qualification  of  this  air¬ 
craft  type.  The  principal  problem  encountered  when  performing  aeroelastic  investigations 
on  aircraft  with  active  control  systems  consists  in  the  fact  that  such  elastomechanic 
structures  cannot  be  regarded  as  (quasi-) linear  and  (quasi-) conservative  systems;  thus 
the  limits  of  validity  of  many  conventional  test  and  analytical  procedures  are  greatly 
exceeded . 

As  a  solution  to  the  various  problems  encountered  during  the  aeroelastic  qualification,  a 
method  for  the  performance  of  ground  vibration  tests  on  aircraft  with  active  control  sys¬ 
tems  has  been  described.  Thereby  the  active  hydraulic  actuators  are  replaced  by  passive 
dummy  actuators.  By  this  step  the  elastomechanic  aircraft  structure  is  linearized  and  re¬ 
veals  conservative  properties  such  that  conventional  ground  vibration  test  methods  can  be 
applied.  Furthermore  an  experimental-analytical  method  for  determination  of  the  nonlinear 
and  preload-dependent  frequency  responses  of  the  actuator  transfer  functions  has  been 
pointed  out  which  permits  the  consideration  of  aircraft  active  control  systems  in  the  gen¬ 
eralized  aeroelastic  equations  by  addition  of  correction  terms. 

The  validity  and  applicability  of  the  elaborated  aeroservoelastic  advances  has  been  checked 
using  the  example  of  a  wing  model  structure  with  a  rudder  driven  by  a  hydraulic  actuator. 

By  comparison  of  the  measured  and  calculated  dynamic  response  of  the  rudder  hinge  angle 
(for  a  given  excitation  of  the  wing) ,  good  agreement  could  oe  shown  between  the  measured 
and  calculated  results.  Furthermore  it  has  been  presented  that  recoupling  effects  due  ei¬ 
ther  to  the  mechanical  systems  (steering  mechanisms)  or  to  the  electric  feedback  system 
(active  control  system)  can  widely  affect  the  dynamic  behavior  of  the  aircraft.  It  is 
therefore  of  primary  importance  to  consider  these  feedback  systems  very  carefully  when 
performing  aeroelastic  qualification  tests  or  analyses. 
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